BACKLUND TRANSFORMATIONS AND DARBOUX INTEGRABILITY 
FOR NONLINEAR WAVE EQUATIONS 



JEANNE N. CLELLAND AND THOMAS A. IVEY 



Abstract. We prove that second-order Monge-Ampere equations for one function of two 
variables are connected to the wave equation by a Backlund transformation if and only if 
they are integrable by the method of Darboux at second order. 



1. Introduction 

Roughly speaking, a Backlund transformation is a method for generating new solutions 
for a given partial differential equation by starting with a 'seed' solution to the same (or a 
different) PDE and solving an auxiliary system of DDEs. Backlund's original example was 
a transformation that produced new pseudospherical surfaces from old, and it is equivalent 
to the following system: 

Vx - Ux = \sm{{u + v)/2), 

^ ' (1.1) 
Vy + Uy = —2 sin((M — v)/2). 

Given an arbitrary smooth function u{x, y), this overdetermined system for v is inconsistent. 
However, if u satisfies the sine-Gordon equation Uxy = sinw then the system is consistent, 
and the function v{x,y), determined up to a constant of integration, will also satisfy the 
sine-Gordon equation. The transformation works in reverse, too: given a solution v{x,y) for 
sine-Gordon, the system determines a 1-parameter family of solutions u{x,y) for the same 
PDE. 

It is this type of Backlund transformation — connecting solutions of two second-order 
Monge-Ampere PDE in the plane, not necessarily the same equation — which is the gen- 
eral subject of this paper. (A second-order Monge-Ampere equation is a PDE where the 
highest-order derivatives may appear nonlinearly but only in the form of the determinant of 
the Hessian.) Another important example of this type is the system 

Zx - Ux = -2exp{{u + z)/2), 

Zy + Uy = exp{{u - z)/2). 

In this example, if z{x,y) satisfies the wave equation (in characteristic coordinates, z^y = 0), 
then the system determines a 1-parameter family of solutions of Liouville 's equation Uxy = 
e^, and conversely. Backlund transformations of this subtype — where one of the two PDE 
involved is the standard wave equation — are the specific concern of this paper. 

Liouville's equation also has the rare property that it is Darboux-integrable — in other 
words, it can be solved by the method of Darboux. (This will be defined below.) The main 
point of this paper is that this is not a coincidence; more precisely, we will prove 
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Theorem 1. Let (M^,X) be a hyperbolic Monge- Ampere system. If there is a normal 
Backlund transformation with 1- dimensional fibers linking this system with the wave equa- 
tion, then the first prolongation of X is Darboux-integrable. Conversely, if the first prolon- 
gation of X is Darboux-integrable, then near any point p G M there is an open set ?7 C M 
around p such that the restriction of X to t/ is linked to the wave equation by a normal 
Backlund transformation. 

The technical terms in this theorem must be explained. Any single PDE or system of PDE 
may be re-cast as an exterior differential system (EDS) or differential ideal (i.e., an ideal, with 
respect to wedge product, in the ring of differential forms on a manifold, that is also closed 
under the exterior derivative), in a way that solutions are in one-to-one correspondence with 
submanifolds to which all the differential forms in the EDS pull back to be zero. (These 
submanifolds, which must usually also satisfy a nondegeneracy condition, are known as 
integral manifolds of the EDS.) In particular, a Monge- Ampere equation in the plane can be 
re-cast as the following type of EDS: 

Definition. A Monge-Ampere exterior differential system is a differential ideal X on a 5- 
dimensional manifold M, such that near any point of M, X is generated algebraically by one 
1-form 6 and two 2-forms Qi,Q2- (Hence, d6 must equal a linear combination of the fi's, 
plus possibly a wedge product with 6 as factor.) The 1-form 6 is required to be a contact 
form, i.e., 6 A d9 A d6 ^ 0. A Monge-Ampere system is hyperbolic if the fi's may be chosen 
so that both are decomposable. 

For example, for Liouville's equation we may take M to be with coordinates x, y, u,p, q, 
and let 

9 = du — pdx — qdy, Qi = {dp — e^dy) A dx, Q2 = {dq — e^dx) A dy. (1.3) 

(Note that d6 = —Qi — Q2-) Given a solution u = f{x,y) of the PDE, we can construct a 
surface S C such that i*6 = 0, i*Qi = i*Q2 = (where z : E is the inclusion map) 

by setting u = f{x,y), p = fx{x,y) and q = fy{x,y). Conversely, any surface S satisfying 
i*6 = 0, i*Qi = i*Q2 = and the nondegeneracy condition i*dx A dy ^ is the graph of a 
solution constructed in this way. 

In the body of the paper, we will also use another type of EDS: 

Definition. A Pfaffian exterior differential system is a differential ideal X on an arbitrary 
manifold M, defined by a vector bundle / C T*M, such that a differential form belongs to 
X if and only if it is a linear combination of wedge products involving sections of / or their 
exterior derivatives. (In practice, our Pfaffian systems will be specified by giving a list of 
1-forms that span the fiber of / at each point.) The rank of a Pfaffian system is the rank of 
the vector bundle. 

A Pfaffian system satisfies the Frohenius condition or is said to be integrable if the exterior 
derivative of any section of I is in the algebraic ideal generated by /. Any Frobenius system is 
locally equivalent to a (possibly underdetermined) system of ordinary differential equations; 
see Chapter 1 in [S]. 

Theorem [T] is about relations between exterior differential systems; in particular, we have 
the following definition from [8]: 
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Definition 1.1. A Bdcklund transformation between two exterior differential systems (M,X) 
and (M,X) is a manifold "B equipped with submersions vr : S ^ M and vf : 23 ^ M (see 
diagram below) and vector bundles J, J C T*!B such that 

(i) the fibers of vr and W are transverse in 

(ii) the rank of J equals the dimension of the fibers of tt, and sections of J pull back to the 
fibers of vr to span the cotangent space of each fiber; 

(iii) J is similarly related to the fibers of vf; 

(iv) the algebraic ideal 3 generated by 7r*T and sections of J is the same as the algebraic 
ideal generated by W*I and sections of J, and is a differential ideal. 



The impact of the last condition is that if !N" C M is an integral manifold of X, then sections 
of J pull back to 7r^^(!N) to satisfy the Frobenius condition, so that integral manifolds of d 
inside 7r~^(!N) may be constructed by solving ODE; moreover, the image under W of each of 
these integral manifolds is an integral manifold of X. Because the definition is symmetric, this 
also works in the other direction: given an integral manifold of X, we can solve a Frobenius 
system on the inverse image in B to obtain a family of integral manifolds of X. For example, 
given a solution z{x,y) of the wave equation, substitution in (11.21) gives an overdetermined 
system of ODE for a solution u{x, y) of Liouville's equation, and in this context the Frobenius 
condition is exactly the compatibility condition for the ODE system. 

The condition of normality for Backlund transformations, assumed in Theorem [H will be 
explained in §2. 

A hyperbolic Monge-Ampere system is a special case of hyperbolic EDS: 

Definition. A hyperbolic EDS of class /c is a differential ideal defined on a manifold of 
dimension k + 4 that, near any point of manifold, is generated algebraically by k 1-forms 
and two decomposable 2-forms. 

Associated to a given hyperbolic EDS X of class k are two characteristic distributions, 
one corresponding to each decomposable 2-form generator. At each point, the distribution 
is given by the 2-dimensional subspace of the tangent space annihilated by the k 1-forms of 
the system and the factors of the chosen decomposable 2-form. (These annihilators form a 
rank k + 2 Pfaffian system, known as a characteristic system of X.) A hyperbolic EDS X is 
integrable by the method of Darboux, or Darboux-integrable for short, if both characteristic 
distributions have two independent first integrals, i.e., functions which are constant along all 
curves tangent to the distribution, and whose differentials are pointwise linearly independent 
from the 1-forms of X. Such functions are also known as characteristic invariants, since they 
are constant along integral curves of the distribution. 

The Darboux-integrability condition has the virtue that it is easy to check, using only 
differentiation and linear algebra, by calculating the successive derived systems of each char- 
acteristic system. An extensive discussion of hyperbolic EDS and Darboux-integrability, 
with worked-out examples, is available in Chapter 6 of [S]. For the purposes of this paper, 
we will need a few more facts about the method of Darboux: 
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• It is known that any Monge- Ampere system which is Darboux-integrable (i.e., has 
two characteristic invariants for each distribution) is equivalent to the standard wave 
equation under a contact transformation (see, e.g., Thm. 2.1 in [2j). 

• If a hyperbolic Monge- Ampere system X has a pair of independent first integrals for 
exactly one of its characteristic distributions, then X is said to be integrable by the 
method of Monge or Monge-integrable for short. (The analogous term for hyperbolic 
EDS of arbitrary class k is Darboux semi-integrable.) 

• If a hyperbolic EDS of class k fails to be Darboux-integrable, it is possible that its 
prolongation, which is a hyperbolic EDS of class /c -|- 2, is Darboux-integrable. Thus, 
a given Monge-Ampere system may lead to a hyperbolic system that is Darboux- 
integrable only after sufficiently many prolongations. 

Prolongation of an EDS is essentially the process of adding higher derivatives as new vari- 
ables and adjoining to the ideal the differential equations satisfied by the higher derivatives. 
For example, for Liouville's equation we add variables r and t to stand for u^x and Uyy re- 
spectively, and adjoin the 1-forms Oi = dp — rdx — e^dy, O2 = dq — e^dx — tdy. The new ideal 
is a Pfaffian system on R'' (with coordinates x,y,u,p,q,r,t) generated by 1-forms 9o,9i,92- 
(Note that the 2-forms Qi and ^2 given in (II. 3p are now in the ideal generated algebraically 
by these 6'o, 6*1, 6*2.) The set of algebraic generators of the new ideal is completed by com- 
puting the exterior derivatives of 61,62, and these are expressible as linear combinations of 
the decomposable forms Q'l = {dr — pe^dy) A dx, Q2 = {dt — qe^dx) A dy, modulo multiples 
of 60,61,62- (Thus, the new ideal is a hyperbohc EDS of class 3.) Let Ai, A2 be the corre- 
sponding characteristic distributions for the prolongation (i.e., Aj is annihilated by 60,61,62 
and the factors of VL[). To see that the system is Darboux-integrable, note that x,r — 
are invariants for Ai and y,t — are invariants for A2. 

Remark. Both Darboux-integrability and the transformation (11.21) enable one to express 
all solutions of Liouville's equation via specifying two arbitrary functions and integrating 
systems of ODE, and these two solution methods are equivalent, although Darboux's method 
requires one to solve more ODEs. For, as mentioned above, substituting the wave equation 
solution z = f{x) + g{y) in (11.21) produces two compatible ODEs for u{x, y). Given an initial 
value for u, these may be integrated simultaneously in the x- and ?/-directions to propagate 
a solution over an open set in the x?/-plane. On the other hand, under Darboux's method 
we obtain ODEs by setting one invariant in each characteristic system to be an arbitrary 
function of the other, yielding in this case the equations 

Px - = Hx), qy - \q^ = i){y), 

which, together with Ux = p and Uy = q, may be integrated to obtain the solution. (In fact, 
the data for these two methods are related by = /" and = g", but in other cases we 
cannot expect the relationship to be this simple.) 

Next, we will put Theorem[T]in context with other results both classical and modern. Much 
of what was known in the 19th century about solving second-order PDE for one function of 
two variables was summarized in Goursat's treatise [B]. In Volume 2, §181 of that work, we 
find the following result: 

Theorem 2 (Darboux-Goursat). Suppose that a second-order PDE for z as a function of 
X, y has the property that there exists a Pfaffian system 

dFi = ^ida + ^idl3, l<i<(., 
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and formulas 



x = V^, y = V2, z = Vs, (1.4) 

where $j, Vi, V2, V3 are functions of Fi, . . . ,Fi, a, /?, /(a), and finitely many of the 
derivatives of / and g, such that (11.41) satisfies the Frobenius condition for arbitrary choices 
of functions / and g, and gives an implicit solution of the PDE for arbitrary choices of initial 
data for the Frobenius system. Then the PDE is Darboux-integrable after finitely many 
prolongations. 

The hypotheses of the Darboux-Goursat theorem are fulfilled if the given PDE is linked 
to the standard wave equation by a Backlund transformation. (For, the d'Alembert formula 
gives solutions of the wave equation Uai3 = in the form u = f{a) + g{j3) for arbitrary / 
and g, and the Pfaffian system in the theorem is given by the equations of the Backlund 
transformation.) Compared with the Darboux-Goursat theorem, one direction of our the- 
orem has a stronger hypothesis (essentially, that £ = 2 and only first derivatives of / and 
g are involved) and a stronger conclusion (that at most one prolongation is required to get 
Darboux-integrability) . 

In Theorem 6.5.14 in |$j it is shown, by an elementary argument, that Darboux-integrability 
of the prolongation implies that there is a Backlund transformation between the prolonga- 
tion (not the original Monge- Ampere system, but one defined by an EDS on a 7-dimensional 
manifold) and the wave equation (defined by an EDS on a 5-dimensional manifold). How- 
ever, this asymmetric transformation-relating the 2-jets of solutions of one PDE to the 1-jets 
of another — is less than satisfying, compared to more symmetrical transformations like (11. 2p . 
Our analysis in §4 shows that it is a much more delicate matter to show that there exists 
a Backlund transformation between two Monge- Ampere systems. We should also note that 
the argument given in [8] for the other direction (Backlund-equi valence to the wave equation 
implies Darboux-integrability) is unfortunately incorrect, and the proof we give in §3 of this 
paper is along completely different lines. 

We now briefly outline the rest of the paper. In §2 we set up the basic machinery re- 
quired for the first half of the proof, namely, the G-structure for Backlund transformations 
originally introduced by the first author in [3]. In §3 we prove the forward direction in our 
theorem by following the implications (for the invariants of the G-structure) of the existence 
of a Backlund transformation to the wave equation. In §4 we prove the converse direction 
by constructing, for any given Darboux-integrable Monge-Ampere equation, an involutive 
exterior differential system whose solutions are such transformations. In §5, we discuss 
our results in the context of earlier classifications of Darboux-integrable equations and of 
Backlund transformations to the wave equation; we also outline an alternate proof technique 
for the converse direction, which can in some examples be used to establish global existence 
of the transformation. In §6 we discuss further steps in our research program. 

We are grateful to the referee who read the first version of this paper, and gave us many 
useful comments and suggestions. 

2. G-STRUCTURE FOR BACKLUND TRANSFORMATIONS 

The material in this section is taken from the first author's paper [3J; additional details 
may be found there. 
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Suppose that (M,X) and (M,X) are hyperbolic Monge- Ampere systems, with 

x = {^,fii,fi2}, x = {^,ni,n2}. 

As a special case of Definition 11.11 we define a Backlund transformation between (M, X) and 
(M, X) to be a 6-dimensional submanifold S C M x M for which the puUbacks to S of the 
forms Qi,Q2,^i, ^2 have the property that 

ni = ni mod{^,^}, i = l,2. 

(The vector bundles J, J G T*'B mentioned in Definition 11.11 are in this case spanned by 
the puUbacks of 6 and 6, respectively.) A Backlund transformation is normal in the sense 
of Theorem [1] if the puUbacks to S of the 2-forms d6, dO are linearly independent modulo 

{ere}- 

Now let ^ be the ideal on B generated by the puUbacks of X and X; according to the 
conditions above, S is generated algebraically by the forms {6, 9, d9, d9}. 

Since X and X are hyperbolic, locally there exist 1-forms , u'^ , , on S such that 
{9, 9, u^, up' , up , u;"^} is a coframing of 3 (i.e., a set of 1-forms that restricts, at each point, 
to be a basis for the cotangent space of S) and 

a = {^, 9, uP NuP, up Nuj^}. 

(It is important to note that 9 and 9 are each separately determined up to a scalar multiple, 
since each determines the contact structure on a 5-manifold.) Any such coframing has the 
property that 

d9 = Aiuj'^ Au"^ + A2UJ^ Au'^ mod {9, 9}, 
d9 = Asu^ Auj'^ + A4 uj^ A mod {6*, 9} 

for some nonvanishing functions Ai, A2, A^, A4. Since d9, d9 are required to be linearly 
independent 2-forms at each point of B, we must have A1A4 — A2A3 7^ 0. 

By rescaling the a;* and adding multiples of 9 and 9 to the tu' if necessary, we can arrange 
that 

d9 = Aiuj^ Auj^ + uj^ Auj^ mod 9, (2.1) 

d9 = uj'^ A u'^ + A2 u]^ A mod 9 

for some nonvanishing functions Ai, A2 on S with A1A2 7^ 1. This coframing is not unique; 
any other such coframing {9, 9, uj^, uP ^ uP ^ cu^} has the form 
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611622 ~ 612621 
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where 611622 — 612621 7^ 0, 011022 — 012021 7^ 0. (The inverse is included for greater ease of 
computation in carrying out the method of equivalence.) A coframing satisfying (12.11) is 
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called adapted, and the group G of matrices of the above form is called the structure group of 
the equivalence problem. (In fact, the most general choice of structure group would include 
a discrete component interchanging the distributions {a;^,a;^} and {uj^ju"^}. However, this 
freedom does not contribute anything crucial to the structure group, and it is easier to work 
with a connected group.) The associated G-structure is the principal G-bundle V 
whose local sections are precisely the adapted coframings over a neighborhood of S. 
In [3], it is shown that V has structure equations 
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(2.3) 



'6 A (A1C2 uj^ - AiGi cj2) + Ai cj^ A + A u^' 
e A {A2C4 - A2C3 u^) + A + ^2 A uj^ 
BiO AO + Ciuj'^ Auj^ 
826 AO + G2Uj'^ Auj^ 
BsO AO + G3UJ^ Auj^ 
B^e AO + G4UJ^ Auj^ 

for some functions Ai,Bi, Gi and 1-forms a,. Pi on V. These equations are chosen so that the 
matrix in (12.31) takes values in the Lie algebra q of G; this is a standard step in the method 
of equivalence. (See [5] for details.) 

The 1-forms ai,Pi are linearly independent from each other and from 6,6,00'^; they are 
called pseudoconnection forms, or more concisely (but imprecisely) connection forms on V. 
They are well-defined only up to transformations of the form 



ai I— > ai + ri uj^ + r2 tu^, 
0:2 0:2 + r2 uj^ + ra , 
0:3 (— > 0:3 + r4 uj^ — ri up' , 



^ — riLj 



r2<jJ 



P\^ P\^ s\ + S2 UJ^, 
P2 ^ P2 + S2 + S3 uj"^, 
Ps^ P3 + Sa uj^ - Si u^, 

Pi ^ Pi — Si — S2 ijj'^. 



(2.4) 



Remark. The coefficients Ai,Bi, Gt are called torsion functions. They may be interpreted 
as the components of well-defined tensors associated to the Backlund transformation, as 
follows. 

A hyperbolic Monge- Ampere system naturally equips the underlying manifold with a 
line bundle L C T*M and two rank 3 characteristic bundles Ki, K2 C T*M whose intersec- 
tion is L. (The generator 1-form ^ is a section of L, and the factors of the decomposable 
generator 2- forms Vti and VL2 span a complement of L within Ki and K2, respectively.) The 
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G-structure for the Backlund transformation shows that is equipped with a well-defined 
splitting of its cotangent bundle: 

T*'B = L®L®Wi®W2, (2.5) 

where L and L are the pullbacks of the Monge- Ampere line bundles from M and M re- 
spectively, and Wi, W2 are spanned by 00^,00"^ and uj^,uj'^ respectively. The normal Backlund 
condition implies that 

7r*Ki = L®Wi, 7r*Ki = L®Wi, i = 1, 2. 

The splitting ( 12. 5p induces a corresponding splitting of A^T*!B, whose summands include 
the line bundles L ^ L, A'^Wi and A'^W2- We may then define a natural map from sections 
of L to sections of A'^Wi, given by 

5o : 9 i-^ projection of d9 into A^H^i. 

But this map is linear under multiplication by functions, and so gives a well-defined map 
between the corresponding vector bundles. The structure equations (12. 3p show that Aiuj^Auj'^ 
is the value of Sq applied to the first member 6 of the coframe. Hence, Ai is the component, 
with respect to the give coframe, of a well-defined tensor in L* ® A'^Wi. Similarly, A2 is the 
component of a well-defined tensor in L ® A'^W2 
We may similarly define a map on sections of Wi by 

61 : u ^ projection of du into Li (g) L2, 

which again is linear under multiplication by functions. The structure equations show that 
Bi9 A 9 and B29 A 9 give the value of 61 on the basis sections u^, u? of W\ respectively. Thus, 
the vector \Bx B-^ gives the components, with respect to the coframe, of a well-defined 
tensor in ®Lx®L2- In a similar way, we see that [i?3 B^ are the components of a tensor 
in W2 ®Lx® L2, [Ci C2] are the components of a tensor in W* A'^W2, and [C3 C4] are 
the components of a tensor in W2 ® A'^Wi, all defined by taking the exterior derivative of a 
section and projecting into the appropriate summand of A^T*B. 

The following results are proved in [3]: 

Proposition 2.1. The vectors [Bi B2], [B3 B4], [Ci C2], [C3 C4] are relative invariants: 
given any point m G B, they are each either zero for every adapted coframing at m, or 
nonzero for every adapted coframing at m. 

Proposition 2.2. If [Ci _C2]_ = [C^ C^] = [0 0], then [Bi B2] = [B3 B4] = [0 0] 
as well, and (M,X) and (M,X) are each contact equivalent to the Monge- Ampere system 
representing the standard wave equation. 

Proposition 2.3. If [Ci C2] = [0 0] (resp^[C3 C4] = [0 0]), then [Bi B2] = [0 0] (resp., 
[B3 B4] = [0 0]) as well, and (M,X) and (M,X) are each Monge-integrable. 

Proposition 2.4. Suppose that the vectors [Ci C2] and [C3 C4] are both nonzero. Then 
the vectors [Bi B2] and [-B3 B^] are either both zero or both nonzero. 

If [Bi B2] = [B3 B4] = [0 0], then the differential ideal generated by {u^ , u"^ , , u"^} 
is a Frobenius system. (The converse is true as well.) It follows that locally, there exists a 
4-manifold V which is a quotient of S and for which the 1-forms u^, u;^, u^, uj^ are semi-basic 
for the projection p : S — V. (Here "locally" refers to the fact that any point in B has a 
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neighborhood which possesses such a quotient, and "semi-basic" means that the restrictions 
of the cij* to the fibers of the projection vanish identically. See [S] for details.) In fact, this 
quotient factors through each of the quotients tt : —>■ M. and vf : S ^ M, as shown by the 
following commutative diagram. 



IT 



TT 



P M 



V 

The vanishing of the vector [Bi B2] implies that the span of {00^,00"^} is unchanged along 
the fibers of p, and is thus the pullback via p of a well-defined rank 2 sub-bundle of T*V. 
When the vector [B^ -B4] also vanishes, the ideal {u^ A cu^, A u"^} is the pullback via p of 
a well-defined hyperbolic system Ti of class on V, and X, X are both integrable extensions 
ofH. 

Backlund transformations of this type are called holonomic. One can test whether a 
Backlund transformation is holonomic by checking whether the Pfaffian system on B spanned 
by the intersection of the basic forms for vr with the basic forms for yf is Frobenius. Note 
that the basic forms for vr are spanned by the Cartan syste of X, and the basic forms for 
W are spanned by the Cartan system of X. 

Holonomic Backlund transformations are generally considered less interesting than non- 
holonomic Backlund transformations because of their limited capacity to generate new solu- 
tions, which we now explain. Given an integral surface K of (M,X), solving the Frobenius 
system 3 on 7r~^(}sf) produces a 1-parameter family of integral surfaces !Na of (M,X). Revers- 
ing the process, beginning with any one of the integral manifolds l^x, in turn produces a 1- 
parameter family of integral surfaces ^A.^t of (M, X). In general, this results in a 2-parameter 
family of integral surfaces of (M,X), and iterating the process produces an ever- increasing 
family of new integral surfaces for both Monge- Ampere systems. 

For example, consider the system fll.2p . If we substitute the trivial solution z{x,y) = 
of the wave equation into fll.2l) . the resulting overdetermined PDE system for u yields the 
1-parameter family of solutions 

u{x,y) = -2\n{-x - ly - ci) (2.6) 

to Liouville's equation. Reversing the process, substituting (12.61) into (11.21] produces a PDE 
system for z which has a 2-parameter family of solutions: 

z{x, y) = 2 ln{-y - C2) - 2 ln(2x + 2ci - C2). (2.7) 

Finally, substituting (12.71) into (II. 2p produces a PDE system for u which has a 3-parameter 
family of solutions: 

u{x, y) = -2 In [c-sxy + (C2C3 - l)x + (ciCi - icaCg - |)?/ + (C1C2C3 - ^cjcs - Ci)) . (2.8) 

It is clear that the solutions (12.61) form a proper subset of the solutions (12.81) . since the 
argument of the latter contains an xy term. 

For the system (II. 2p . and for non- holonomic Backlund transformations in general, succes- 
sive iterations of this process continue to produce new solutions. However, if the Backlund 



The Cartan system for a given EDS 2 is the smaUest Frobenius system that contains 2. 
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transformation is holonomic, tlien all integral surfaces of (M, X) and (M, X) produced by this 
process must lie in the inverse image of a single integral surface of (V, 7i). It follows that 
successive iterations can produce no more than a 1-parameter family of integral surfaces for 
each Monge- Ampere system. 

3. Proof that Backlund implies Darboux 

Now suppose that we have a Backlund transformation as in §21 and that the Monge- 
Ampere system (M^,X) is contact equivalent to the standard wave equation Zxy — 

0. We 

can choose local coordinates (X, Y, Z, P, Q) on M such that X is generated by the forms 

e = dZ - PdX -QdY, ni = dX A dP, ^2 = dV A dQ. 
There is a unique local section cr = (6*, 9, u^, uP', , uj^) : S ^ P satisfying 

B = dZ-PdX-QdY, 
00^ = dX + CiO, 

uj^ = dP + C2e, (3.1) 

00^ = dY + 030, 

u}^ = dQ + de 

for some functions Cj on B. (Note that, because specifying this portion of the coframing 
determines a unique local section of P, the 1-form 9 is also uniquely determined.) The 
functions Cj are the pullbacks under a of the corresponding torsion functions, and this 
coframing has Ai = 1. When the structure equations (12. 3p are pulled back to !B via a, the 
1-forms ai,(3i — which were linearly independent from the 1-forms 6,6,00"^ on V — must pull 
back to some linear combinations of these 1-forms. 

Now we embark on the process of comparing the structure equations (12.31) to those for the 
explicit coframing above. First, note that 

= d{dX) = d{uj^ - CiO) = -{ai A cu^ + as A cu^ + Ciuj'^ A 00^) mod 9. 

Therefore, by choosing r2,r3 appropriately in (12.41) . we may assume that 

ai = aio6 + aiiiu^ + ^Cilu"^, 02 = a2o9 - \CiUJ^ 

for some functions aio,aii,a2o on B. Similar considerations of d{dP),d{dY),d{dQ) modulo 
6 yield similar expressions for the remaining connection forms: 

tti = aiod + aiitu^ + (3i = bioO + bi3Uj^ + ^Csu^, 

012 = a2oO - P2 = hoO - ^C^u^, 

as = asoO + \C2uP' , /^s = &30^ + \C4^'^, 

oiA^ = a^oO — ICscu^ + 0421-1-'^, Pi = b^oO — ^C^uj^ + b^^tu'^. 

Next, a straightforward computation shows that for the coframing (13. ip . 

de = eA {C2UJ^ - Ciuj^ + C4J' - Cacj^) + cu^ A + A a;^ 
Comparing this with the first structure equation in ( 12. 3p yields 

^13 = fC*4; ^44 = ~|C*3- 
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In order to continue this comparison, we need to introduce derivatives of the functions A2, Ci. 
So, set 

dA2 = A2fie + ^2,0^ + A2,lUJ^ + A2,2UJ^ + A2-3UJ^ + A2,4Uj\ 

dC2 = C2,o9 + C2fle + C2,lUJ^ + C2,2CU' + Ca.sCu' + C2,4CU^ (3.2) 

rfCa = C^flO + + Cs,!^^ + C3,2^' + C3,3a;3 + C3,4a;^ 

Comparing the structure equations (12.31) for du'^ with the derivatives of the exphcit forms c<j* 
in (13.11) yields 

Oio = Ci,i — C1C2, biQ = Cs^3 — C3C4 

020 = C*l,2 + C*! ; ^20 = ^3,4 + Cf; 

0-30 = C*2,l — C2, &30 = (^4,3 ~ Cl, 

040 = C*2,2 + C1C2, &40 = ^4,4 + C^Ci, 

in addition to the following relations among the torsion and its derivatives: 



(3.3) 



While we don't have an explicit coordinate representation for 6, we can still explore the 
consequences of d{d9) = 0. Computing d{d9) = modulo 9 yields 

_ -2^2,1 + ^2(^2 + 2) _ -2^2,2 - Ci(A2 + 2) 

2A2 ' 2A2 ' 

A2,0 = ^2(C3,3 + (^4,4 — Cl,l — (^2,2)- 

Then computing d{d9) = modulo 9,00^,00"^ yields 

Ci C2 

^2,0 ~ ^1,2 + ^^2,2 — ^2(^^3,3 + ^^4,4) ^^2,1 ;P^2,2 — ^3^2,3 — ^4^2,4- 

A2 A2 

(Note that the fact that 6' is a contact form implies that A2 cannot be zero.) 

At this point, all coefficients in the structure equations (12. 3p have been expressed in 
terms of the functions y42, Ci, C2, C3, C4 and their first derivatives. In addition, we have 
relations among the derivatives that amount to an overdetermined PDE system for these 
five functions on B. Necessary compatibility conditions for this system may be found by 
computing d{dA2) = d{dCi) = 0. In particular, computing 

d{dCi) = d{dC2) = mod 9, 9, uj\uj^, 

d{dC3) = d{dCi) = mod 9, 9, cu^ cu^ 
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Bi = 




Cl,3 


= C1C4, 


^3,1 


= C2C3, 


B2 = 


~^2,0; 


Cl,4 


= — C1C3, 


^3,2 


= — C1C3, 


B3 = 


~^3,0' 


^2,3 


= C2C4, 


^4,1 


= C2C4, 


B,= 


~^4,0' 


^2,4 


= — C2C3, 


^4,2 


= — C1C4. 



yields 



Ci,o = ~^2Cifi — ^2^1,2 + + Cs.s + C'4,4) 

^2,0 = ~^2C'2,0 ~ C*iC2,l + C2(— 6*2,2 + 6*3,3 + 6*4,4) 

6*3,0 = ^639 — 6*46*3,4 + 6*3(6*1,1 + 62,2 — 6*3,3)) 
6*4,0 = ~6*4o — 6*36*4,3 + 6*4(6*1,1 + 6*2,2 — 6*4,4). 



(3.4) 



At this point, we have derived all the relations among the torsion functions on S and their 
derivatives that will be needed in order to prove that (M, X) is Darboux-integrable after at 
most one prolongation. The proof of Darboux-integrability is divided into two main cases. In 
§3.H we prove Darboux-integrability under the assumption that the vectors [6*1 6*2], [6*3 6*4] 
are both nonzero. This case is further divided into three subcases, depending on the ranks of 
certain Frobenius systems that arise during the proof. Precise statements of the results are 
contained in Propositions 13.11 13.21 and 13.31 In §3.21 we prove Darboux-integrability under 
the assumption that exactly one of the vectors [6*1 6*2] , [6*3 6*4] vanishes; the precise result 
is contained in Proposition 13.41 As noted in ^ it is not necessary to consider the case where 
[6*1 6*2], [6*3 6*4] both vanish, since in that case both Monge- Ampere systems are contact 
equivalent to the standard wave equation. 



3.1. Case 1: [6*1 6*2], [6*3 6*4] are both nonzero. Without loss of generality, we may 
assume that 6*2 and 6*4 are nonzero. Consider the exterior derivatives of the ratios ^ and 
A straightforward computation shows that 



C2 



d (^^^ = mod 9,9,uj\u^; 



therefore, d must lie in the last derived system of JCi = {9, 9,uj^,lj^} — i.e., the largest 
integrable subsystem of /Ci, denoted by }C^^\ Similarly, 



mod 9,9, 00^,00^^, 



so d { ^) must lie in the last derived system of IC2 



{9,9,uj^uj''}. 



First, consider the system /Ci. In order to compute its first derived system )C^i \ we must 
find those 1-forms in /Ci whose exterior derivatives are zero modulo the linear span of the 
1-forms in /Ci. To this end, we compute: 



d9 



d9 = A2UJ^ A uj^ 
du^ = 6*1 A 
du^ = Co 00^ A 00^ 



> mod /Ci. 



Therefore, /Cj^^ = {9- ^2^, u;^ - Ci9, - 6*2^}. Observe that 

cji - Ci9 = dX, uj^ - C29 = dP. 
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Therefore, the rank 2 subsystem {u^ — Ci6, uJ^ — = {dX, dP} of IC^i^ is integrable, 

and the next derived system /c[^^ (i.e., the first derived system of IC\^^) contains this rank 
2 system. The only question is whether, in fact, ICf'^ = /C^^"* — i.e., whether K,^^^ itself is 
integrable. In either case, we will have JCf'^ = }C^^\ A computation shows that 

d(9 - A^e) =eA [(^2,3 + A2CiiA2 - 1)) cu^ + (^2,4 - ^2^3(^2 - 1)) cu^] mod /CS'\ (3.5) 

so the rank of /c[°°'* is either 3 or 2, depending on whether or not the 1-form in brackets 
vanishes. 

Similarly, we can compute that 

/c^^) = (0-e,u}^- c^e, - de} = {e-e, dY, dQ}. 

So K,'^^ contains the rank 2 subsystem {u^ — C^O, uj^^ — C^O} = {dY, dQ}, and a computation 
shows that 



d{9 -9) = -9 A 



(^2,1+^2(^2-1)) , (A2,2-Ci(A2-l)) 2 

UJ H uj 

Ao Ao 



mod/C^^\ (3.6) 



So the rank of IC''^^ is either 3 or 2, depending on whether or not the 1-form in brackets 
vanishes. 

Now we must divide into cases depending on the ranks of these derived systems. 

3.1.1. Case 1.1: /c[°°'* and K,^^ both have rank 3. In this case, we have the following result: 

Proposition 3.1. If [Ci C2], [C3 C4] are both nonzero and /C[°°'* and /Cg""^ both have rank 
3, then the system (M,X) is contact equivalent to the standard wave equation. 

Proof. By Theorem 2.1 of ^, it suffices to show that (M,X) is Darboux-integrable, i.e., that 
each of the characteristic systems {9,00^,00"^} and {9,u!^,u!^} — which are well-defined on M 
even though the 1-forms cu* are not — contains a rank 2 integrable subsystem. 

The hypothesis that /C[°°'' and IC'^^ both have rank 3 implies that the expressions 03.51) 
and (13.61) must both vanish identically; therefore, 

A2,l = -^2(^2-1), ^2,2 = CM2-I), ^2,3 = -A2CM2-I), A2,4 = ^26-3(^-1). 

Using these conditions, a straightforward computation shows that 

(9, U0\ CU^jd) = {A2UJ^ - Ci'9, A2UJ^ - C2^}, 

{9, u\ u^Y^^ = {^207^ - C^9, A2U^ - Ci9}, 
and that each of these derived systems is integrable. □ 

3.1.2. Case 1.2: Exactly one of k!'^^ and ]C^2°^ has rank 3. Without loss of generality, we 
may assume that JC^^ has rank 2 and is equal to {dY,dQ}, and that JC^^^ has rank 3. It 
follows that (13.51) vanishes identically and that (13. 6p does not. Since all our results are local, 
we will assume that we are working on an open set where (13. 6p is nonzero. The vanishing of 
(13. 5p implies that 

A2,3 = -^2^4(^2 - 1), ^2,4 = ^2C^3(^2 " 1)- 
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Recall that the function 7^ satisfies 



d(^^^e}Ci'-^ = {dY,dQ}. 
It follows that ^ is a function of Y and Q alone. Now consider the 1-form 

This 1-form is contained in the span of uj'^,uj'^, and we have 

duj^ = mod uj^; 

so uj^ is a multiple of an exact 1-form, say = XdY. Moreover, because uj^ is expressed 
solely in terms of Y and Q, A and Y may be chosen to be functions depending only on Y 
and Q, and which are therefore well-defined on M. The crucial point here is that there exists 
an exact 1-form in the span of {u^,u'^} which is well-defined on M. Then we have 

dY A dQ = A dQ = XdY A dQ = dY A dQ, 

where 

rQ ^ 

Q{Y,Q)= / \{Y,t)dt. 
Jo 

Since 

de = dX AdP + dY A dQ, 
Pfaff's Theorem (see Ch. 1 of [8]) implies that there exists a function Z on M such that 

e = dZ -PdX-QdY. 

We can now repeat all our constructions starting with the coordinate system (X, F, Z, P, Q), 
but now our adapted coframing a will have the additional property that tu^ = dY and C3 = 0. 
Thus we will drop the tildes and assume that C3 = for the remainder of this subsection. 

Proposition 3.2. If [Ci C2], [C3 C4] are both nonzero, /C[°°'' has rank 3, and K,^^'^ has rank 
2, then the system (M,X) is Monge-integrable, and it becomes Darboux-integrable after one 
prolongation. 

Proof. The same argument as that given in Case 1.1 shows that the characteristic sys- 
tem {^, tu^, cu^} on M contains a rank 2 integrable subsystem; therefore, (M,X) is Monge- 
integrable. 

In order to prove the second half of the Proposition, we will need to make use of relations 
among the second derivatives of the functions Ci, C2, C4. These will be denoted as, e.g.. 

Note that, although (for example) the A2^ij are second derivatives of A2, because we are 
working in a coframing rather than in coordinates, we cannot assume that ^2,ij = ^2,jj- 
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Computing d{dA2) = mod {9, cv^, cv^} shows that 

^2,10 = ^2(C'4,41 — Cl,ll — C2,21 + C*2(C'l,l + (^2,2 — (^4,4)) 

+ ^2, 1(^^4,4 — C2,2 — C1C2) + ^2,2(C'2,1 — C2), 

^2,20 = ^2(C4,42 ~ Cl,12 ~ (^2,22 + C'l(C'4,4 ^ C*!,! ^ C'2,2)) 

+ ^,l(Cl,2 + CD + A,2(C4,4 - + C1C2), 

^2,21 = ^2,12 + (^2 — 1) ( C*!,! + 6*2,2 T^^l r^C2 

V ^2 ^2 

Next, computing d{dC^ = mod 9} shows that 

6*4,34 = 6*4,43 + (^2 ~ 1)64,0 + 64(61,1 + 62,2), 



64,31 — 62(64,3 + 64;, 



2^ 

2\ 



64,32 = —61(64,3 + 64), 

64.41 = 6264,4, 

64.42 = —6164,4. 

Now computing d{dCi) = mod 9 shows that 

r -rr , ((^2 + 1)61-^2,1) ^ 

•-^4,01 — ^4:^2,0 "I '-"4,0' 

^ _ rr ((-42 + 1)^2 + ^,2) ^ 

•-"4,02 — — '-"4<^1,0 1 '-^4,0) 

^2 

^4,30 = ^4,03 " ^26464^0) 
^4,40 — ^4,04) 

and then computing d{dCi) = mod 9 shows that 

64,30 = 64(61,13 + 62,23 — 64,43 + 640) + (64,3 — 6|)(6i,i + 62,2) — 6403, 
64,40 = 64(61,14 + 62,24 — 64,44) + 64,4(61,1 + 62,2) — 64^04, 

61,12 = —62,22 + (^2 — l)6i_o ~ 6164,4 + — ^^4,o> 

62,21 = —61,11 — (^2 — 1)62 + 6264,4 — -^^-^ j-^ ^64 0- 

/I2O4 

Next, computing d{dCi) = mod {9, 9} yields 

61.13 = 64(61,1 + 6162), 

61.14 = 0, 

61.23 = 64(61,2 — 61), 

61.24 = 0, 

61,21 — —62,22 — 26164,4 + 2(^2 — 1)61 + — -7—^ ^64 0- 

/i204 
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Similarly, computing d{dC2) = mod {6, 6} yields 

^2,13 = Ca{C4:2,1 + C^)-, 
^2,14 = 0, 

^2,23 = ^'4(^2,2 ~ CiC2)i 
6*2,24 = 0, 

6*2,12 — "6*1,11 + 2C2C4,4 — 2(^2 — 1)6*2 ~ ~ r~7^ ^^4 0- 

Computing d{dCi) = mod {uj'^ , , u^} yields 

6*1,10 = 6*26*2,22 ~ 6*161,11 — ^26]^ 01 + ((2 ~ ^2)6*2 — ^2,1)6*1,0 + 6*4,4(6*1,1 + 26*16*2) 
((A -1)61-^,2) 



^26*4 



-6*26*40, 



n - I (^2,1 - 6*2) ^ 

^^1,10 — ^1,01 A ^1,0' 



and computing ^(^62) = mod {uj^,uj^,uj^} yields 

6*2,20 — 6*16*1,11 ~ 626*2,22 — ^262,02 + ((^2 ~ 2)6*1 ~ ^2,2)6*2,0 + C*4,4(6*2,2 ~ 26*16*2) 

((A2- 1)62 + ^,1) 



+ 



^26*4 



-6*16*40, 



n — I ("^2,2 + 6*1) ^ 

'-^2,20 — '-^2,02 "T ^ '-^2,0) 



Now computing d{dA2) = mod cu^ yields 



^2,10 



-7- [—6*1^42,11 — 6*2^2,12 + ^2,1(6*16*2 + 6*2,2) + ^2,2( — 6*2,1 + 6*2)] 
^2 

- {A2 - 1)^2,0 - 6*2(6*1,1 + 6*2,2) - {{A2 - 2)6*2 + 2/12,1)6*4,4 



((^2- 1)6*2 + ^2,1) 
^26*4 



6*, 



4,0' 



A. 



2,20 



1 



—C1A2A2 — 6*2^2,22 ~ ^2,1 6*1,2 + 



6*1 



62 



+ ^2,2 6*1,1 - 



6*16*2 



+ (A2 - 1)6*1,0 + /(C^i,i + C2,2) + ((^2 - 2)6*1 - 2^2,2)6*4,4 

, ((^- 1)6*1-^2,2) ^ 



and 



2,14 



A 



2,24 — '-'4,44 



6*4 



0. 



Finally, we need two additional relations which do not become apparent until we differen- 
tiate the equations for dC^^^ and ^6*4,4. Computing (i((i6*4,3) = mod {9 — 9, 00^,00"^ — 6*4^} 
yields 

(6*46*4,03 - ((^2 + 1)6*4^ + 6*4,3)6*4,0)^ A [(^2,1 + 6*2(^2 - 1)) + (^2,2 - 6*1(^2 - 1)) ^']. 
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Note that the right-hand factor is precisely (13. 6p . which we have assumed is nonzero. There- 
fore, 

Q,o(C4.3 + {A2 + l)Cl) 

^4,03 — 

Precisely the same argument applied to d{dCiiA) shows that 

_ _ C'4,oC4,4 

^4,04 — • 

With these relations in hand, consider the characteristic system /C = {6,uj^,uj'^} of X — 
which is well-defined on M, even though uj^ is not. We need to show that after one pro- 
longation, the corresponding characteristic system /C of the prolongation contains a rank 2 
integrable subsystem. In order to perform this computation, we need to construct a basis for 
K, consisting of 1-forms which are well-defined on M. Fortunately, 9 and uj'^ are already well- 
defined on M. For the remaining 1-form, it will be convenient to choose a 1-form which is 

contained in the first derived system /C*'^'' = {a;^ — C4 6*, u^}. To this end, introduce functions 
r, (7 on !B such that the 1-form 

is well-defined on M. (The fact that K^^^ is well-defined on M guarantees the existence of 
such functions.) As before, we denote the derivatives of these functions by 

dr = TqO + Tq6 + Tiiu^ + T2UJ^ + Tsuj^ + t^uj'^, 

dg = goO + go9 + giuj^ + 5-2^^^ + gsuj^ + g^^'^, 

and similarly for second derivatives. 

Because ip is well-defined on M, dtp contains no terms involving 6. This, in turn, determines 
the partial derivatives To,go: 

Tq = 6*4^4, 

90 = Cl — 6*4^3 — 5'C4^4. 

We will also need to make use of relations among the second derivatives of r, g. These are 
determined by computing d^dr) = d{dg) = 0; this is a straightforward computation, which 
we omit here for the sake of brevity. 

We can define a partial prolongation I of X on M x M as follows. (Note that Darboux- 
integrability of the partial prolongation implies Darboux-integrability of the full prolonga- 
tion.) Let t be a new coordinate on the M factor; then the partial prolongation if is generated 
by the 1-forms 6,6 = ip — tuj^, and the 2-form tu^ A u"^. Again, this system is well-defined on 
M X M, even though uj^,uj'^ are not. 

A straightforward computation shows that 

dd' = -TTi A a;^ mod {6, &'}, 

where 

Hi = dt — (tTi — gi) uy" — {tr^ — e'^5'2) '^'^ ■ 
The corresponding characteristic system of x' is 

k1 = {6,6',ni,uj^}. 
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We will now compute the derived systems of /C and show that the second derived system 
/C is a Frobenius system of rank 2; this will complete the proof of the Proposition. In 
order to compute the first derived system, we compute: 



> mod {6',6' ,7ri,u;^} 



de =0 

dixi = Euj^ A 
duj^ = Q 

(the last line following from C3 = 0), where 

E = e^Cl + (e"^4 - tu)C^ - e^C^,^ - {e^ g + t)C4,4 + {e^ % - tr^). 
Let = TTi — E6, so that 



Next we compute the derived system of JC 
—I C4 — 



/(I) 



d9 



A. 



-9 A [(^2,1 + ^2(^2 - 1)) + (^2,2 - - 1)) 



F ^ 



d'K2 = A [(^2,1 + ^2(^2 - 1)) + (^2,2 - Ci{A2 - 1)) U^] 



duj^ 
where 







mod {9 , 712, o;^}, 



Once again, we see the bracketed 1-form in fl3.6p appearing. Since this 1-form is assumed to 
be nonzero, the derived system /C has rank 2 and is spanned by the forms 00^ and 

VTs = C47r2 - F9. 
Finally, another computation shows that 

diT^ = 0" 



duj^ = 



mod {7T3,uj }; 



therefore, /C is integrable. 



□ 



3.1.3. Case 1.3: /C[°°'* and IC'^^ both have rank 2. Now we assume that the bracketed 1-forms 
in both (13. 5p and (13.61) are nonzero. By the same argument as that given in the previous 
case, we may assume that Ci = C3 = 0, with = dX, = dY. 

Proposition 3.3. If [Ci C2], [C3 C4] are both nonzero and IC^^^ and IC'^^ both have rank 
2, then the system (M,X) becomes Darboux-integrable after one prolongation. 
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Proof. The proof is very similar to that of Proposition I3.2[ We must now consider both 
characteristic systems 

ICi = {e, uj\ u^}, IC2 = (e, u^} 

oil. As before, these systems are both well-defined on M, even though cu^ and are not. 
We introduce functions p, r, f,g on "B such that the 1-forms 

7] = eP{A2Uj'^ - 026 - fuj^), ip = e^{uj^ - - gu^) 

are well-defined on M. These forms have the property that 

We construct the prolongation x' of X on M x as follows. Let r, t be new coordinates 
on the factor; then the prolongation if is generated by the 1-forms 9,6i = rj — ruj^, 62 = 
ip — tu^, and their exterior derivatives. 

The remainder of the proof consists of applying the argument of Proposition 13.21 to each 
of the characteristic systems /C'^, IC'2 of the prolongation if . The argument varies only in the 
details of the calculations, and so we omit it for the sake of brevity. □ 

3.2. Case 2: One of the C-vectors vanishes. Without loss of generality, assume that 
[Ci C2] = [0 0], and that d ^ 0. By Proposition ESI it follows that [Bi B2] = [0 0] as 
well. 

Proposition 3.4. If [Ci C2] = [0 0], then the system (M,X) is Monge-integrable, and 
it becomes Darboux-integrable after at most one prolongation. Furthermore, the Backlund 
transformation 3 C M x M is holonomic. 

Proof. It follows from Proposition 12.31 that (M,X) is Monge-integrable; in fact, the charac- 
teristic system {9,uj^,u!^} contains {uj^,u^} = {dX,dP} as a rank 2 integrable subsystem. 

Now consider the other characteristic system JC = {9,u^,u'^} of X. One easily computes 
that the first derived system of JC is 

In order to find the second derived system, we compute: 

d{uj^ - CsO) = ^eA {A2,lUj' + A2,2UJ^) ] 

P_ \ mod^W. 

d{uj^ - C^9) = -^eA iA2,lCu' + A2,2UJ^) 

A2 ) 

If A2^i = ^2,2 = 0, then /C*^^-* is integrable; in this case, X is Darboux-integrable and hence 
contact equivalent to the standard wave equation. Therefore, we assume that ^42,1 and ^2.2 
are not both zero. 

In order to prove the second statement, we will construct a partial prolongation of X and 
proceed as in §3.1.2[ But first we need to derive relations among the derivatives of the torsion 
functions. 

Picking up where we left off at (13.41) . consider d{dCs),d{dC4). Computing d{dC^) = 
d{dCi) = mod {e -e,uj^- C^e, - C^} yields 

^3,0 = C*4o = 0. 
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From fl3.3p . it follows that 3^ = 3^^ = 0; therefore, the Backlund transformation is holonomic, 
as claimed. 
We now have 

dc^ = c,,^{uj^ - a,e) + C3,4(^' - c,e) = c^^.dv + c,,,dQ, 

dC^ = C4,3(^^ - CsO) + C4,4{uj^ - C^e) = C^^sdY + d^^dQ. (3.7) 

It follows that C3, C4 are functions of Y and Q alone. Now the same argument as that given 
in §3.1.21 shows that we may assume C3 = 0; moreover, C4 remains a function of Y and Q 
alone when we do so. Computing d{dCi) = yields the following relations among the second 
derivatives of C4: 

'-'4,30 — — L'4<-'4,34; L^4,40 — —0404^44, 04^43 — 04^34, 

^4,30 = C'4,40 — C'4,31 = ^4,32 = ^^4,41 = ^4,42 = 0. 

Now consider the characteristic system /C = {^, cj^, cu'*}. As we computed above (recalling 
that C3 = 0), its first derived system is 

/cw = {cu^cu^-C4^}. 

As in §3.1.21 choose functions (7, r so that the 1-form 

is well-defined on M, and construct the partial prolongation x' of X and the 1-form Q as we 
did there. Similar calculations to those of §3.1.21 show that the corresponding characteristic 
system /C' of X' has a rank 2 integrable subsystem. This completes the proof. □ 



4. Proof that Darboux implies Backlund 

4.1. The non-Monge-integrable case. In this subsection (M,X) is assumed to be a hy- 
perbolic Monge- Ampere system which is Darboux-integrable after one prolongation, but not 
Monge-integrable. We will construct a canonical coframing associated to the prolongation. 
We will then use this coframing to construct an integrable extension ^ of (M, X) in such a way 
that % defines a Backlund transformation between (M,X) and the standard wave equation 
ZxY = 0. 

Lemma 4.1. Near any point of M, there exists a coframing (6*, tti, 712, r/^, r/^) such that 6 
spans the 1-forms of X, and the characteristic systems Ci,C2 of X have derived flags 

Proof. By a result of Juras |10j , (M, X) is locally contact equivalent to a system encoding a 
PDE of the form 

u^y = F{x,y,u,p,q). 

Thus, there are local coordinates x, y, u, p, q near the given point of M such that X is generated 
by the 1-form 6 = du — pdx — qdy and the 2-forms {dp — F dy) A dx and {dq — F dx) A dy. 
It is easy to verify that the coframing given by 6, = dx, rf = dy, tti = dp — F dy — Fq6, 
and 712 = dq — F dx — Fp6 has the properties claimed. □ 
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In terms of the local coframing on M given by the lemma, the prolongation (M',X') is 
defined as follows: let M' = M x M^, with coordinates r, t on the factor, and let I' be the 
Pfaffian system on M' generated by 6 and the forms 

e^='jri-rr]\ Q^ = n2-trf. (4.1) 

Lemma 4.2. Near any point of M' there exists a coframing (^, 6*1, 6*2, r/""^, r/^, vts, 1x4) such that 
X' is generated by 9, 01,02, satisfying 

d0 = -01 A r]^ - 6*2 A mod 

d0^ = -TTa A 7]^ mod 0, 0i (4.2) 
d02 = —TT4 A mod 0, 02, 

with the derived flags of the characteristic systems of X' given by 

C[ = {0,0,,02,V\^3} 3 {0,0i,v\^3} ^ {0i,v\^3} 3 {^,^3} = C[^^\ 
C'2 = {0,0i,02,ri\'r^A] 3 {^,^2,r/^7^4} D {^^2,r7^7^4} D {7/^,714} = C'2^^\ 

Proof. Let 6*, 77^ = (ix, r^^ = dy be part of the local coframing on M (pulled back to M') 
constructed in the proof of Lemma WA\. and let 0i, 02 be defined as in (14.11) . Then 

d0i = -{dr - {DxF)dy) A dx mod 0, 0i, 
d02 = -{dt- {DyF)dx) A dy mod 0, 02, 

where 

D^F = F, + Fup + Fpr + F^F, 
DyF = Fy + F^q + FpF + F^t. 

For the moment, let 7T3 = dr — {DxF)dy. Because dn^ = modulo dx, tc^, 0, 0i, 02, it follows 
that TTs lies in C'i^^\ Moreover, we may subtract a multiple of from tts to ensure that vts 
lies in C^^^^ 

Next, we prove that the last derived system of C[ has rank 2, rather than rank 3. (A 
similar argument applies to C!^.) Suppose that C^*-^-* = {0i,r]^,7i^} is integrable. From (14. ip . 
it is clear that this is equivalent to the statement that {711,77^,713} is integrable — i.e., that 
dTTi = dr]^ = dn^ = modulo ni, rj^, tc^. But tti and rj^ are both well-defined on M, so their 
exterior derivatives do not involve 7r3. It follows that diTi = dr]^ = modulo 7ri,?7^, and 
C^^ = {7ri,?7^} is integrable, contrary to the hypothesis that (M,X) is not Monge-integrable. 

□ 
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The conditions in Lemma 14.21 are preserved by changes of coframing of the form 
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with ai, a2, c ^ 0. Let G C GL{7, M) be the group of such transformations, and let V be the 
G-structure on M' of which the coframing of Lemma 14.21 is a section. 

After absorbing as much torsion as possible and differentiating to uncover relations among 
the torsion, V has structure equations 
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2C202 A TXi 





(4.4) 

Because of the dimensions of the derived flags of the characteristic systems (given in 
Lemma [4.2p . Ai,Bi are not both zero, and ^2,-62 are not both zero. Furthermore, we can 
choose a local section a : M' — P satisfying the conditions that rj^ = dx, rf = dy, and the 
forms vr3,7r4 are integrable; i.e., 

dn^ = mod 773, d7[4 = mod tc^. 

To see why, note that {77^,773} is a Frobenius system, and so it is spanned locally by two 
exact 1-forms. Thus we can adjust tt^ by adding multiples of rj^ in order to make it a multiple 
of an exact form. Similarly, we can add multiples of 77^ to tt^ in order to make a multiple 
of an exact form. However, we cannot independently scale 713 and 774 to make both of them 
exact. 
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This choice of section is not unique; it is determined up to a transformation of the form 
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with c 7^ 0. However, we can make the choice of a unique (albeit shghtly non-canonical) as 
follows: since vts, 714 are integrable 1-forms, we must have 

TTa = e^fi^i, 7r4 = e'^d^2 

for some functions ^1,^2,1,9 on M'. Using the remaining scaling freedom, we can arrange 
that h = —g; the resulting coframing a : P — *• M is uniquely determined. 

When we pull back the structure equations via a, the pseudoconnection forms ai, a2, Pi, 
f32, 7 become semi-basic. By making use of the remaining ambiguity in these forms and the 
conditions imposed thus far on the coframing, we can assume that 

ai = {D^ + E^)r]' 

"2 = (^2 + ^2)^^ 
Pi = 2D^ns 
P2 = 2D2vr4 

7 = -CiO^ - C2Q2 - Elv' - E^v^ + i^lTTg + F27r4 
for some functions Ci, Di, Ei, Fi. Then the structure equations for this coframing become: 

dd = dh {-Cidi - C2Q2 - Elf]' - E2f + FiTTg + F27Ti) - 9i A T]^ - 02 A f 
dOi =01 A (-0202 - E2Tf + FiTTg + ^3774) + Di0i A 7]^ - A T]^ + A {A202 + B2rf ) 

d02 = 02 A {-Ci0i - Eiri^ + FiTTg + F2TXi) + 0^02 A rf - Tx i A rf + A (Ai^i + BiT]^) 
dr]^ = (4.6) 
df = 

dns = 713 A {Ci0i - C202 + EiT]^ - E27f + ^3774) 

dn^ = 71^ A {-Ci0i + C202 - EiT]^ + E2Tf + FiTTs). 

(Note that these torsion functions are completely unrelated to those in ^and ^) 

As in ^ we will need to compute relations among the derivatives of the torsion functions 
in order to show that (M,X) has a Backlund transformation to the wave equation. We begin 
by differentiating the structure equations (14.61) . Using notation similar to that in §31 we 
denote derivatives as, e.g., 

dAi = Ai^oO + Ai,i0i + Ai,202 + Ai,3V^ + ^iaV^ + ^i.svrg + Ai^ir^. 
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(Note that since this coframing is defined on a different manifold from that in ^ the indexing 
of the derivatives is different as well.) 

Computing d{d6) = d{d6i) = d{d92) = d{dTTs) = d{dTT4) = yields the following equations 
for the derivatives of the torsion functions: 

dA^ = Ai,o^ + ^1,1^1 - 3AiC202 + ^1,3^' - Ai{D2 + 2E2W + + AiF27r4 

dA2 = ^2,0^ - SAaCi^i + ^2,2^2 - ^2(^1 + 2^1 )r7^ + ^2,4^/^ + ^2^^371-3 + ^2^2714 
dBi = Bi^oO + (Ai,3 - AiDi)ei - 251^2^2 + Bi^^r]^ - Bi{D2 + ^2)//' + A^^ 
dB2 = B2fl6 - 252^1^1 + (^2,4 - ^2^2)^2 - B2{Di + Ei)r/^ + ^2,4//' + ^2714 

dCi = AiC2e + Ci^iOi - 010202 + Ci,3r/^ - + O1E2W + Ci,57r3 + CiF27r4 

d02 = A2O19 - Oi029i + C2,2^2 - (1^2 + C2^i)r/^ + 02,iV^ + 02Fins + C2,6vr4 (4.7) 
dD, = B,02e + + §^2^2 + Di^sv' + 1(3^2 - 5i)r/2 + (2Ci - ^1,5)713 

dD2 = B2O19 + lAiOi + ^2,2^2 + i(35i - ^2)??' + ^2,4^?' + (2C2 - ^2,6)7r4 

= B,02e + (Ci,3 - OiDi)e, + \A2e2 + ^i,3r/' + \{B2 - Bi)r]^ + E,,,7r^ 
dE2 = B2O19 + lAiOi + (^2,4 - C2 1^2)^2 + \{Bi- B2W + E2,aV^ + ^2,6Vr4 

dFi = (2CiFi - Ci,5)^l - C2Fi^2 + (Cl + ElFi - ^1,5)77! - ^2^17^2 + Fi,57r3 + Fi,67r4 
dF2 = -^1^201 + (2C2F2 - C2,6)^2 - EiF2r]^ + {O2 + ^2i^2 - ^2,6)^' + F^fiT^Z + i^2,6Vr4. 

Because Ai appears as a derivative of Bi, and cannot vanish simultaneously, Bi 

cannot vanish on any open set in M'. In fact, Bi cannot vanish identically on any fiber of 
the projection M' M, and the same is true of B2. Henceforth we restrict to the dense 
open set in M' where B2 are both nonzero, and note that this set surjects onto M. 

We may obtain further relations among the derivatives of the torsion functions by differ- 
entiating equations (14. 7p . Computing d{dAi) = d{dBi) = modulo 6*, 6*1, r^^, 714 yields 

^1,0 = ^1(^2,4 ~ -^2,2 ~ O2D2), Bi Q = -81(6*2,4 — -D2,2 — O2D2). 

Then computing d{dBi) = modulo 6,61,1]^ yields C2,6 = O2F2. Similar considerations of 
d{dA2) and d{dB2) show that 

^2,0 = ^2(Ci,3 — -Di,i — CiDi), 52,0 = -82(^1,3 — -Di^i — Ci-Di), 6*1,5 = CiFi. 

It will now be convenient to derive several equations and solve them simultaneously. First, 
d{dOi) = modulo 9i,ri^ implies that 

v4i(C2,2 - 0^2) = MCi,i - O!) (4.8) 
Ai(3C2,4 - D2,2 - SO2D2) = 252(Ci,i - Ol). (4.9) 

Additionally, d{d02) = modulo 62,7]^ implies that 

^2(3^^1,3 - - 3CiDi) = 25i(C2,2 - O^). (4.10) 

Finally, d{dDi) = modulo 6*1,77^,713 implies that 

fil (3^2,4 - D2,2 - 3O2D2) = 52(3Ci,3 - /^i,i - 3CiDi). (4.11) 
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The general solution to equations (I4.8l) - (l4.1ip is most easily expressed in terms of a new 
torsion function if, such that 

= - A^H, = 3(Ci,3 - CiZ^i) + 2B^H, 

C2,2 = C*2 ~" ^^H, 1^2,2 = 3(C2,4 ~ C2D2) + 2B2H. 

Next, we need information about the derivatives of Ci^3 and 6*2,4. Computing d{dCi) = 
modulo 9i yields 

rfCi,3 = (^2^1,3 - iAiA2 + B^C,C2)9 + Ci,3l^^l + (IA2C1 - ^2^1,3)^^2 
+ C,,33V' - 1(^1,3 + 2E2C,,3 + - 52))r/2 

+ (FiCi,3 - Ci^i,5 + 2C! - A^H)7rs + F2C^,^7i^, 
and computing d{dC2) = modulo 62 yields 

C^C'2,4 = {CiA2^4 — ^AiA2 + B2CiC2)0 + (IA1C2 — CiC2,4)^l + C2,42^2 
- i(/l2,4 + 2^1^2,4 + C72(52 - fii))r/l + ^2,44^/' 

+ i^lC2,47r3 + (F2C2,4 - ^^3^2,4 + 2Cl - A2H)Tl4. 

Now, computing d{dCi) = d{dC2) = 0, d{dDi) = modulo ri^,TC3, and d{dD2) = 
modulo r/^, 714 yields four different expressions for dH. Taking linear combinations of these 
expressions shows that 

A(A,3 - A^D{) = BMi,i - AlCi) (4.12) 
^2(^2,4 - A2D2) = ^2(^2,2 - A2C2). (4.13) 

Equations fl4.12p and (14.131) may be solved by introducing new torsion functions Ji, J2, such 
that 

= AiCi + AiJi A2,2 = A2C2 + A2J2 

Ai^3 = AiDi + BiJi A2,A = A2D2 + B2J2. 

Then the various expressions for dH may be combined to show that 

Ci,3i = 2CiCi,3 - C^,D, + A,C2 - IA1J2 - {AiDi + BiJi)H 

^2,42 = 2(726*2,4 ~ C2D2 + A2C1 — 1^2 Jl — (A2D2 + B2J2)H, 

and 

dH = {HiCi^s - CiD^ + B^H + C*2,4 - ^72^2 + B2H) - \C1J2 + C2J1) 6 

+ CiHOi + C2He2 + {DiH + \J2W + (^2i^ + \Ji)rf + Fiif7r3 + ^21/714. 
They also imply the relation 

i/(Ci,3 - C*iDi + 5ii7) - \C1J2 = H{C2,4 - C2D2 + B2H) - \C2J1. (4.14) 
The equations for dAi, dA2 now take the form: 

dAi = -2Ai{C2,4 - C2D2 + B2H)d + Ai{Ci + Ji)di - 3AiC2^2 

+ {AiDi + fii JOr]' - Ai(D2 + 2^2)7/2 + AiF37r3 + ^1^2774 (4.15) 
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dA2 = -2A2(Ci,3 - C^D^ + B^H)9 - SAaCi^^i + ^2(^2 + J2)^2 

- A2iDi + 2Ei)r]^ + {A2D2 + ^2 J2)^' + ^2i^3vr3 + A2F2TT4. (4.16) 
Now, computing d{dBi) = modulo 61, yields 

^2,44 = (-D2 — E2)C2,A + C2D2A — HB2^A — {D2 + E2)B2H + C2D2E2 + B2{Ci — Ji). 

Similarly, computing d{dB2) = modulo 6*2,77^ yields 

Cl,33 = (^1 - ^l)Ci,3 + CiDi,3 - i/fil,3 - {D, + + CiDi^i + 5i(C2 - J2). 

Next, computing d{dAi) = shows that 

dJi = 4A^C29 + {CiJi + AiKi)9i- C2J162 

+ (2Ci,3 - + 25i/f + 2^2,4 - 2^2 ^2 + 2^2// + Ji + fiii^Or^i 

- (2Ai + E2Ji)rf + Fi Ji7r3 + F2 Ji7r4 
for some function /^i. Similarly, computing d{dA2) = shows that 

dJ2 = 4A2C19 - CM + (C2 J2 + A2K2)92 - (2^2 + E,J2)v^ 

+ (2Ci,3 - 2CiDi + 2BiH + 2^2,4 - 2C2D2 + 2B2H + D2J2 + 82X2)1^ 

+ F1J27T3 + F2J2TT4 

for some function i^r2. But now computing d{dH) = modulo 9, 9i, 92 yields 

B^{K^+4H) = B2iK2 + 4H). 

It follows that 

Ki = -AH + B2M 
K2 = -AH + BiM 

for some function M. Computing d{dH) = modulo 92, rj^ and d{dH) = modulo 6'i,?7^ 
shows that 

AiM = A2M = 0. (4.17) 

Claim. M = 0. 

Proof. Suppose not. Then by fl4.17p . Ai = A2 = 0. Therefore, equations fl4.15p and (14.161) 
reduce to 

= dAi = BiJiT]^ 

(} = dA2 = B2J2V^. 

Since Bi, B2 are nonzero, it follows that Ji = J2 = 0. Then 

= dJi = (2Ci,3 - 2CiDi - 2BiH + 2C2A - 2C2D2 + 2^2^ + BiB2M)r]^ 
= dJ2 = (2^1 3 - 2CiDi + 2BiH + 2^2,4 - 2C2D2 - 2B2H + BiB2M)Tf. 

Subtracting the two coefficients above yields 

4(^2 - 5i)i7 = 0, 

so either Bi = B2 ot H = 0. 

First suppose that i?i = i?2- Then 

= d{B2 - Bi) = -2Bi{Ci9i - C292) mod ^,r/\r/2, 
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so Ci = 6*2 = 0. It follows that 6*1^3 = 6*2,4 = as well. But now 

= dJi = B^M, 

so M = 0, as desired. 

Now suppose that H = 0. Computing d{dCi^3) — Cid{dDi) = modulo 9 yields 

BIB2M = 0, 

so M = in this case as well. □ 
Finally, computing d{dH) = 0, keeping the relation (14.141) in mind, yields two additional 
relations: 

2C2(2Ci 3 - + 2B^H) + (2C2 - J2)(2C2 4 - 2C2D2 + 2B2H) + ^2(4^1 - Ji) = 

(4.18) 

(2Ci - Ji)(2Ci,3 - 2C1D1 + 2B^H) + 2Ci(2C2,4 - 2C2D2 + 2B2H) + A^^AC^ - J2) = 0. 

(4.19) 

We now have all the relations that will be needed for the involutivity calculation below. 

Now suppose that M = carries a Monge-Ampere system I representing the wave 
equation Zxy = 0, generated algebraically by the contact form 

9 = dZ-PdX-QdY (4.20) 

and the 2-forms dP AdX and dQAdY. If there were a Backlund transformation 3 C M x M, 
then Z would be a local coordinate on the fibers of 3 — > M and the functions X, Y, P, Q on 
S would satisfy the Backlund condition 

{dPAdX,dQAdY} = {7riAri\7i2Av^} mod 9,9 (4.21) 

(see the definition at the beginning of §2). 

Accordingly, we let S = M x M, with coordinate Z on the second factor. We will show 
that, on an open neighborhood of any point of "B, there exist functions X, Y, P, Q such that 
the ideal ^ = TU {9} on B (where 9 is defined as in (14.201) ) gives a Backlund transformation 
between (M,X) and (M,X). We will do this by setting up an EDS whose integral manifolds 
correspond to functions satisfying these conditions; once we know that this EDS is involutive, 
an application of the Cartan-Kahler Theorem will prove the existence of the desired Backlund 
transformations. 

Let B' = M' X M, again with Z as the coordinate on the second factor; we extend the 
projection M' ^ M to a projection "B' —>■ hj the identity on the second factor. It will 
be convenient to set up our EDS in terms of the coframing {9, 9, 9i, 92, f]^,"!]"^, 7^3, 7^4) on B'. 
Thus, we will regard X, Y, P, Q as functions on B', but require that 

dX, dY, dP, dQ e {9, 9, 9i, 92, 7]\ r]^} 

so that they are in fact well-defined on B. In order to satisfy the Backlund condition (I4.2ip . 
we will furthermore require that 

{dX,dP}c{9,9,9i,r]^}, {dY,dQ} C {9,9,92,r]^}. (4.22) 

From this, and the structure equations ( 14. 6p . it follows that {dX,dP} (resp., {dY,dQ}) is 
the largest integrable subsystem of {6', 6*, 6*1, 77^} (resp., {9,9,92,r]'^})- Therefore, 

= dx e {dX, dP}, rf = dy e [dY, dQ}, 
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and by a contact transformation on M, we may assume that X = x, Y = y. Thus, we will 
set 

9 = dZ — P dx — Q dy, 

and condition (I4.22p becomes 

dPe {9,9,61,7]'}, dQ e{9,9,92,rf}. 

Suppose that 

dP = Po9 + P^9 + Pi9i + PsT]' (4.23) 
dQ = Qo9 + Qq9 + Q292 + Q^v^. (4.24) 

Observe that normality of the Backlund transformation requires that Pi,Q2 ^ ^ and Pi ^ 

Remark. Equations f l4.23p -f H:.24l) give an overdetermined system of first-order partial differ- 
ential equations for functions P and Q. The process of generating compatibility conditions 
for such systems can be carried out systematically by computing the exterior derivatives of 
the 1-form equations, and using the fact that the repeated exterior derivative of a function 
is zero. Moreover, applying Cartan's Test for involutivity (see [8], Chapter 7) to the result- 
ing EDS will tell us when we can stop differentiating: if the system is involutive then no 
further compatibility conditions arise through differentiation, and solutions exist that may 
be constructed by applying the Cartan-Kahler Theorem. 

Differentiating (I4.23p modulo 9,9,9i,ri' yields 

(Po + P0Q2) 92^r]^ = 0, 
and differentiating (I4.24p modulo 9, 9, 92, rf yields 

(Qo + QoPi) OiAv' = 0. 

Therefore, because the 1-forms 9,9,9i,92,i]',i]'^ are linearly independent on S, we have 
Po = — P0Q2, Qo = —QqPi, and we may write 

dP = Pq(9 - Q29) + Pi^i + Psr^i (4.25) 
dQ = Qo{9 - Pi9) + Q2^2 + QaV^. 

Note that neither Pq nor Q^ can vanish identically: for, if Pg = 0, then differentiating 
f l4.25p shows that Pi = as well, which contradicts the hypothesis of normality. (A similar 
argument applies to Qq.) 
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Differentiating equations fl4.25p modulo various combinations of 1-forms leads to the fol- 
lowing expressions for the exterior derivatives of Pq, Pi, P3, Qq, Q2, Qa- 

dP^ = -PoQ-o^e + Po-o(9 - Q^e) + p-,,e, + P-,,r,' 
dP, = P,,i9 - Q,9) - iMl^^^e^ _ c,P,9, 



dPs = {P^, + P§){9-Q2t 



Qo 

, {B,Q2 + QoiiD^ + E{)P^ + 2Po(g2 - ^i))) , 



+ PssT]^ + PlTTs (4.26) 

dQo = -QoPmd + Qmio - PiO) + Qo2^2 + Qmv' 
dQ. = Q-UO - m - C,Q29, - {^^^Ell^MAe, 

dQ. ^ (Qo. + Q§)(^ - P.9) - iB.P.^hm^E2)Q2 + 2Q,iP,-Q2))) ^^ 

Pq 

+ Q^av'^ + (527r4. 

Now computing d{dPi) = modulo 9 — Q29 and d{dQ2) = modulo 6' — Pi9 yields 

^00 = q-p2q^Iq, - p^) (2^1(^1 - ^Q2)PiQ2P^Qo + 2C2PIQ2P0QI 

+ 2(Ci,3 - C^i^i + B^H)P^Q2PoQo + A^{P, - Q2)QlP^ + ^2Pfg§) 

^50 = p p n\ (^C,P,QlP^Qo + 2C2(g2 - 2P,)P,Q2PoQl 

+ 2(^2,4 - C2D2 + B2H)P,QlP^Q^ + Aig3p_2 + ^2(^2 - Pi)PiQt) 
_ Po{A,Q2 + 2C,QoP^) 

g^A ^^^^ 



02 



-^03 



04 



go(A2Pi + 2C2Pog2) 

P^2 

Po{BiQ2 + {Di + E^)P^Q^ + (2g2 - Pi)P^o) 
PiQo 

Qo{B2Pi + {D2 + E2)Q2Po + (2Pi - g2)Pogo) 
Q2P0 



This leaves only P33 and g44 as undetermined second derivatives of P and Q. 

We are now ready to set up our exterior differential system. Let S = 3 x ]R^°, with 
coordinates P, Q, Pq, Pi, P3, Qq, Q2, Qa, P33, g44 on the M}^ factor. Let W be the rank 8 
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Pfaffian EDS on "B generated by the 1-forms 

e, = dp- p^(e - Q^e) - p^e^ - p^ri\ 

Q2 = dQ- Qo(9 - P^e) - Q292 - QaV^ 

and similar forms 63, . . . , 63 prescribing conditions fl4.26p for dpQ, dPi, dP^, dQ-Q, dQ2, and 
dQ4^, substituting the values (14.271) for Pqq, Pqi, Pq^, Qqq, Qq2, Qm- Integral manifolds of W 
satisfying the independence condition 9 A 9 A 9i A 62 A r]^ A r]"^ A A tt^ are in one-to-one 
correspondence with the desired functions P, Q defining a Backlund transformation. 
The structure equations for this EDS have the form: 
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(4.28) 



where 



and \E'i, . . . , are 2-forms which are quadratic in the forms 9, 9, 9i, 92, 1]^, rf, tts, 7r4, with 
coefficients which are polynomial functions of the two quantities 

{2HPi{Q2 - Pi)Q^+Q2{2CiQ2 + Ji{Pi - Q2))Qo + A,Ql) P^ 



- 2C2PIQI + 2(C2,4 - C2D2 + B2H)Qj^P,Q2, 

{2HQ2{Pi - Q2)P^ + Pi{2C2Pi + J2{Q2 - Pi))Po + A2P^)Qi 

- 2CiQlP^ + 2(Ci,3 - C,D, + B,H)Pj^P,Q2. 



(4.29) 



(4.30) 



In order to find integral manifolds, we must restrict W to the locus Z G defined by 
the simultaneous vanishing of (I4.29P and (I4.30p . (These relations could also be found more 
directly, by computing d{dPi) = d{dQ2) = 0.) While relations of this sort could easily lead 
to hopeless incompatibility for the PDE system given by (14.250 and (I4.26p . it turns out that 
differentiating these quantities yields no new relations. 

A case-by-case analysis, based on the vanishing or non-vanishing of various torsion coef- 
ficients of (M,X), shows that the functions Pg, Qq, Pi, Q2 are all nonzero on an open subset 
G Z which is surjective for the projection 2. — > !B. Since normality requires that these 
functions be generically nonzero, we further restrict W to this open subset. 

Let W denote the pullback of W to W is a rank 6 Pfaffian EDS on Z^. Because 
differentiating the equations defining Z yields no new relations, W is torsion-free; moreover, 
it is straightforward to check that W is involutive with last nonzero Cartan character si = 
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2. (See |8] for a discussion of Cartan's test and involutivity.) Therefore, it follows from 
the Cartan-Kahler Theorem that local integral manifolds exist and are parametrized by 2 
functions of one variable. 

We summarize this result as: 

Proposition 4.3. Let (M,X) be a hyperbolic Monge- Ampere system which is not Monge- 
integrable, and which is Darboux-integrable after one prolongation. Then near any point 
p G M there is an open set U C M around p such that the restriction of X to f/ is linked 
to the wave equation by a normal Backlund transformation; moreover, the set of all such 
Backlund transformations is parametrized by 2 functions of one variable. Up to contact 
transformations, all such Backlund transformations preserve the space of independent vari- 
ables x,y. 

4.2. The Monge-integrable case. In this subsection (M,X) is assumed to be a hyperbolic 
Monge- Ampere system which is Monge-integrable and Darboux-integrable after one prolon- 
gation. As explained below, we will construct a canonical coframing associated to the partial 
prolongation of X, and then proceed as in §4.1[ 

A similar argument to that of Lemma 14.11 can be used to prove: 

Lemma 4.4. Near any point of M, there exists a coframing {6, ir i, 7r2, rj^ ,ri'^) such that 6 
spans the 1-forms of X, and the characteristic systems Ci,C2 of X have derived flags 
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Indeed, the same coframing as that given in the proof of Lemma l4?T] satisfies the conditions 
of Lemma [4. 4[ Note that this lemma only assumes the Monge-integrability of (M,X). 

In terms of the local coframing on M given by the lemma, the partial prolongation (M',X') 
is defined as follows: let M' = M x M, with coordinate r on the M factor, and let X' be the 
Pfaffian system on M' generated by ^, the 1-form 61 = 711— rrj^, and the 2-form Ti2 /\ff- 

Lemma 4.5. Near any point of M' there exists a coframing {6,6i,7T2,r]^,i]'^,TTs) such that 
X' is generated by 9, 9i, and tt2 Ai]'^, satisfying 

d0 = -6*1 A r?^ - 6*2 A 77^ mod 6 

(4.31) 

dOi = -Ti^ /\rj mod^,^i, 
with the derived flags of the characteristic systems of X' given by 

C'2 = {9,e,,7r2,v'} 3 {e,n2,v'} 3 {7^2, v'} = C'2^^^- 

Proof. As in the proof of Lemma 14.21 the usual or "full" prolongation of X on M X R2 is 
generated by 6*, 61 and O2 = 1^2 — iff ■ We may construct the coframing (6*, 6*1, 6*2, f]^,f]'^, tts, H4) 
on M X precisely as in Lemma 14.21 and this coframing satisfies the structure equations 

The hypothesis of Darboux-integrability implies that the characteristic system 

}c, = {e,e,,e2,vi,7r3} 

of the prolongation contains a rank 2 Frobenius system. As in Lemma 14.21 it follows from 
the structure equations and the construction of tts that 

IC^^^ = {9,9^,r],,ns}. 
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However, this system is well-defined on M' = M x M, and it coincides with the characteristic 
system C[ given in the statement of the present Lemma. It follows that we may adjust tt^ so 
that it lies in the rank 2 Frobenius system C[^'^\ 

Note that the second characteristic system C2 is simply the sum of X' and the pullback of 
the characteristic system C2 of X, and the structure of its derived flag follows from that of 
C2. □ 

The conditions in Lemma 14.51 are preserved by changes of coframing of the form 
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(4.32) 



with ai, 02, c 7^ 0. Let G C GL{6, M) be the group of such transformations, and let V be the 
G-structure on M' of which the coframing of Lemma 14.51 is a section. 

After absorbing as much torsion as possible and differentiating to uncover relations among 
the torsion, V has structure equations 
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(4.33) 

Because of the dimensions of the derived flags of the characteristic systems (given in Lemma 
14. 5p . A,B are not both zero. Furthermore, we can choose a local section a : M' — > V 
satisfying the conditions that t]^ = dx, yf = dy, tt2 is exact, and ir^ is integrable. The 
resulting coframing is uniquely determined. 

When we pull back the structure equations via a, the pseudoconnection forms ai, 02, 
P21 Ps, 7 become semi-basic. By making use of the remaining ambiguity in these forms and 
the conditions imposed thus far on the coframing, we can assume that 



ai = Eri\ a2 = 0, /?i = Gn2, P2 = 0, ps = {-H+2E)t:^, 7 = FT:2+Gif 
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for some functions E, F, G, H. Then the structure equations for this coframing become: 

de = eA (Ftts + Gr]^) - 01 A T]^ - n2 A T]"^ 
dOi = 9iA (Ftts + Er]' + Gt]^) - tts A r^^ + A {An^ + Br]^) 
d7i2 = (4.34) 
dr]^ = 
= 

d7l3 = 713 A {-C9i + F7l2 + Hf]' + Gt]^). 

Once again, we will need to compute relations among the derivatives of the torsion func- 
tions in order to show that (M, X) has a Backlund transformation to the wave equation. We 
begin by differentiating the structure equations fl4.34p . We denote derivatives as, e.g., 

dA = AqO + AiOi + ^2712 + ^3//^ + A^rf + ^5713. 

Computing d{d9) = d{d6i) = d{d7[^) = yields the following equations for the derivatives 
of the torsion functions: 

dA = Ao9 + ACOi + A27r2 - AEr]^ + A^rf 
dB = BoO + BGOi + A^Tii - BEt]^ + B^t]"^ 
dG = GiOi + CF7r2 + Gsv' + CGt]^ + 6*5713 

dE = EiOi + 2A7r2 + E^t]^ + 2BTf - Gtx^ (4.35) 

dF = AGO + ^2712 + At]^ + F^rf 

dG = BGe + F47r2 + Bt]^ + G^r]"^ 

dH = {GE - G-i)9i + ^712 + H-^T]^ + BTf + H^n-^. 

We may obtain further relations among the derivatives of the torsion functions by differ- 
entiating equations fl4.35p . Computing d{dF) = d{dG) = modulo 7i2,v'^ recalling that 
A, B cannot vanish simultaneously yields 

Ao = A{G,-GE), Bo = BiGs-GE), Gi = -G\ G^ = Q. 

Then computing d{dA) = d{dB) = modulo 6, 7:2,1]'^ yields 

E, = 2iGE-G,), 

and d{dG) = implies that 

dG3 = (G^E - 2GGs)ei + {AG + FG3)n2 + Gs^r]^ + (BG + ^^3)77^ - G^n^. 
Finally, computing d{dA) = d{dB) = modulo 7r2,?7^ yields 

C33 = GE3 + EG3. 

We now have all the relations that will be needed for the involutivity calculation below. 

Now suppose that M = carries a Monge-Ampere system X representing the wave 
equation Zxy = 0, generated algebraically by the contact form 

9 = dZ-PdX-QdY (4.36) 
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and the 2-forms dP A dX and dQ A dY . As in §4.11 if there were a Backlund transformation 

5 C M X M, then Z would be a local coordinate on the fibers of S ^ M and the functions 
X, y, P, Q on S would satisfy the Backlund condition 

{dP ^dX,dQ ^dY} = {'Kl^r]^,■K2^ri^} mode.e. (4.37) 

As in §4.H let B' = M' x R, again with Z as the coordinate on the second factor; we 
extend the projection M' M to a projection B' — B by the identity on the second factor. 
We will regard X, y, P, Q as functions on B', but require that 

dX, dY, dP, dQ G {e,e, Oi, 712, v\ 

so that they are in fact well-defined on B. In order to satisfy the Backlund condition fl4.37p . 
we will furthermore require that 

{dX,dP} C {6,6, 01,7]^}, {dY,dQ} C {6,6,7T2,r]^}. (4.38) 

The same argument as that given in §4.11 shows that by a contact transformation on M, 
we may assume that X = x. However, the same is not true for Y: the system {6,6,TT2,ri^} 
on B contains a rank 3 integrable subsystem, so we cannot necessarily arrange to have 
1]'^ G {dY,dQ}. There are three different, geometrically natural conditions that we could 
impose on the intersection of the rank 2 PfafRan systems {dY,dQ} and {tt2,ti'^}, each of 
them potentially leading to a different type of Backlund transformation: 

(1) {dY,dQ} n {7r2,?7^} has rank 1 and is spanned by a non-integrable 1-form. 

(2) {dY,dQ} n {7^2, v"^} has rank 1 and is spanned by an integrable 1-form. 

(3) {dY, dQ} n {772, V } has rank 2. 

In cases (2) and (3) we can arrange that Y = y via contact transformations on M and M, 
but in case (1) this is not possible. 

4.2.1. Case (1). In this case we have 

6 = dZ - Pdx-QdY, 

and condition fl4.38p becomes 

dP G {6, 6, 6^, ri^}, dY, dQ G {6, 6, 712, r/'}. 

Suppose that 

dP = P06 + Pq6 + P161 + P.r]' 

dQ = Q06 + Q^6 + Qsvra + Q^V^ (4.39) 
dY = Y06 + Y^6 + Y2712 + Y^r]\ 
Normality of the Backlund transformation requires that Pi 7^ 0, Q2Y4 — Q4Y2 7^ 0, and 

Pi + Q2Y^ - QaY2. 

The argument proceeds in much the same fashion as that of §4.11 differentiating equations 
fl4.39p leads to relations among the derivatives of P, Q, Y . Eventually we are led to a Pfaffian 
exterior differential system W whose integral manifolds satisfying the independence condition 

6 A 6 A 61 A 7!-2 A T]^ A T]"^ 7^ are in one-to-one correspondence with the desired functions 
P, Q, Y defining a Backlund transformation. This EDS is involutive with last nonzero Cartan 
character S3 = 1. Therefore, local integral manifolds exist and are parametrized by 1 function 
of three variables. 
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If we impose the additional condition that the Backlund transformation be holonomic, 
we find that the resulting EDS is involutive with last nonzero Cartan character S2 = 2. 
Therefore, among the Backlund transformations of this type, there is a small, proper subset, 
parametrized by 2 functions of two variables, consisting of holonomic transformations. 

4.2.2. Case (2). In this case, we can use contact transformations on M and M to arrange 
that {dY, dQ} fl {712, yf} is spanned by rf = dy = dY. Then we have 

6 = dZ — P dx — Q dy, 

and condition fl4.38p becomes 

dP e{9,9,9^,ri^}, dQ e {9,9,n2,v^}. 

Suppose that 

dP = Po9 + P^9 + Pi9i + P^T]', (4.40) 

dQ = Qo9 + Qq9 + QaTTs + QiV^. 

Normality of the Backlund transformation requires that Pi, (^2 7^ and P\^ Qi- 

Differentiating equations fl4.40p leads to relations among the derivatives of P, Q, and to a 
Pfaffian exterior differential system W whose integral manifolds satisfying the independence 
condition 9 f\ 9 f\ 9i f\ 112 f\ t]^ f\ if 7^ are in one-to-one correspondence with the desired 
functions P, Q defining a Backlund transformation. This EDS is involutive with last nonzero 
Cartan character S2 = 1. Therefore, local integral manifolds exist and are parametrized by 
1 function of two variables. 

If we impose the additional condition that the Backlund transformation be holonomic, 
we find that the resulting EDS is involutive with last nonzero Cartan character si = 3. 
Therefore, among the Backlund transformations of this type, there is a small, proper subset, 
parametrized by 3 functions of one variable, consisting of holonomic transformations. 

4.2.3. Case (3). In this case, we can use contact transformations on M and M to arrange 
that 7f = dy = dY, tt2 = dQ. Then we have 

9 = dZ — P dx — Q dy, 

and condition (14.38P becomes 

dP G {9,9,9i,r]^}. 

Suppose that 

dP = Po9 + Pq9 + Pi9i + P-i7]\ (4.41) 

Normality of the Backlund transformation requires that Pi 7^ and Pi 7^ 1. 

Differentiating equation f l4.4ip leads to relations among the derivatives of P, and to a 
Pfaffian exterior differential system W whose integral manifolds satisfying the independence 
condition 9 A 9 A 9i A tt2 A t]^ A r]"^ 7^ are in one-to-one correspondence with the desired 
functions P defining a Backlund transformation. 

The involutivity calculation in this case depends on the torsion functions in the structure 
equations fOD . If 

AG - BF = AC + F{CE - C3) = AF^ + FA2 - AF2 = BF^ + FA^ - AF^ = 0, (4.42) 

then W is involutive with last nonzero Cartan character S2 = 1, and so local integral man- 
ifolds exist and are parametrized by 1 function of two variables. Otherwise, there are no 
solutions with Pi 7^ 0, and hence no normal Backlund transformations of this type. 
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Observe that in this case, the G-structure on the Backlund transformation S (cf. Q will 
satisfy the condition that (omitting obvious pullback notations) 



{u;',iu'} = {n2,v'} = {dY,dQ}. 



Therefore, all transformations of this type satisfy the hypotheses of Proposition 13.41 and so 
are holonomic. 

We summarize these results as: 

Proposition 4.6. Let (M,X) be a hyperbolic Monge- Ampere system which is Monge- 
integrable, and Darboux-integrable after one prolongation. Then there exist Backlund trans- 
formations of types (1) and (2) above between (M,X) and the standard wave equation (M,X), 
and of type (3) if the torsion functions of (M,X) satisfy (14.421) . The generic Backlund trans- 
formation is of type (1) and does not preserve the space of independent variables. There are 
both holonomic and non-holonomic Backlund transformations of types (1) and (2), and all 
Backlund transformations of type (3) are holonomic. 



In this section we review the classifications of second-order Darboux-integrable Monge- 
Ampere equations, due to Goursat and Vessiot, and discuss the connection between our 
results and the work of Zvyagin. We will also give examples of a method for explicitly 
solving for Backlund transformations linking these equations to the wave equation. 

5.1. The Goursat- Vessiot List. Goursat [7] studied non-linear PDE of the form 



which are Darboux-integrable at the 2-jet level, classifying them up to complex contact 
transformations that preserve the form fl5.ll) . Using Lie-theoretic techniques, Vessiot [llj 
reproduced Goursat 's classification, expanded to include linear equations, and showed that 
some of the equations on Goursat 's list were equivalent under more general contact transfor- 
mations. Recently, Biesecker [1] re-proved Vessiot 's classification using Cartan's method of 
equivalence, with respect to real contact transformations. Retaining Goursat's numbering, 
the list is: 



5. Examples 



xy 



F{x,y,u,u^,Uy) 



(5.1) 



(x + y)u, 



'xy 




(I) 




(III) 

(IV) 



(11) 



u u 



'xy 



±(l){u^)ll){Uy), 



where (j){t),ilj{t) satisfy the ODE df/dt ±t/ f = K for some nonzero constant K; 



{x + y)u^y = 7(M^)7(My), 
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(V) 



where 7 is implicitly defined by 7(t) — 1 = exp(t — 7(t)); 



Ur — u 



ay 


/ 11 \ 


(VI) 




= eVl + 


(VII) 


yuxy 


= /(X, Uxy), 


(VIII) 


Uxy 




(IX) 


^xy 




(X) 


^xy 


= + UxUy. 

\u + X u + y J 


(XI) 



(In (VI) and (VIII) the function / is arbitrary.) To Goursat's original list, Vessiot added 
representatives of the two equivalence classes of Darboux-integrable linear equations: 



u 



xy 



a{x, y)ux + &(x, y)uy - a{x, y)b{x, y)u, (XII) 



where h{x,y) = —a^ and k{x,y) = —by must satisfy the system {\iah)xy = 2h — k, {In k)xy 
2k — h with h k] and finally, 



2u 

{x + yy 

In the above list, we have replaced Goursat's original versions of (VII) and (XI) by simpler 
equations that Vessiot showed were equivalent to them by contact transformations; see [11], 
part 2, pages 5 and 6, respectively. Vessiot also observed that (VI) is contact-equivalent to 
(X), and (VIII) may be reduced by a contact transformation to the special case 

Uxy = (VIII*) 

X + y 



Our Theorem [H together with Goursat's classification, implies the following 

Corollary 3. If a second-order Monge- Ampere PDE for one function of two variables 
is linked to the standard wave equation by a normal Backlund transformation with 1- 
dimensional fibers, then the PDE is either equivalent to the wave equation by a contact 
transformation, or equivalent to one of the equations (I)- (XIII) in the above list. 

5.2. Zvyagin's List. Zvyagin [14j investigated second-order Monge- Ampere equations linked 
to the standard wave equation by a Backlund transformation, and asserted that all such 
transformations that are non-holonomic are exhausted by a list of six examples in addition 
to Liouville's equation. Zvyagin did not publish a proof of this classification, and did not 
give explicit forms for the Monge- Ampere equations for some of the transformations on his 
list. He did give an explicit transformation for Goursat- Vessiot equation (I): 

Vp-v^=^/^, v^-v^=-i/^, (Z.I) 

\j x-y \l x-y 

where, as in §4, Z is the solution to the wave equation, with x- and y- derivatives P and Q. 
(The X- and y-coordinates are preserved by the transformation.) 
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Corollary [3] implies that every one of Zvyagin's transformations must be identifiable with 
an equation on the Goursat-Vessiot list. We have calculated explicit forms for certain trans- 
formations on Zvyagin's list, and we can identify the following transformations as belonging 
to equations (II), (III), and (VII), respectively: 



zp-up = Vz^ - mVi + ZQ-uq = Vz"^ - + <5^; (z.ii) 

(sinh(;z) + |e-^)P + e^^-'^^/VP^ - 1 (sinh(z) - ie-"')Q + e-H+^)/V<5^ - 1 



P = : , Q 



— sm M — sm u 



(Z.III) 



where w is related to u and Z by cosm = coshZ — |e 



p={l- 2e"+^)P - 2e("+^)/Ve"+^ - l/P^Tl, q = -Q - e(«-^)/Ve"+^ - 1- 

(Z.VII) 

Each of the above transformations, which preserve the x- and y-coordinates, may be verified 
as being non-holonomic. To see how this is done, suppose that the transformation equations, 
when solved for p and q, take the form 

P = fix,y,u,Z,P), q = g{x,y,u,Z,Q). (5.2) 

The Cartan system for the Monge- Ampere equation is spanned by dx, dy, du, dp and dq, 
while the Cartan system for the wave equation is spanned by dx, dy, dZ, dP, dQ. Recall from 
§2 that a transformation is holonomic if the intersection of these systems is Frobenius. In 
light of the transformation equations 05.21) . the intersection of these two systems is spanned 
by dx, dy, dp — fu du and dq — Qu du. The last two 1-forms are congruent modulo dx and dy 
to 

^1 = fpdP + fzdZ, ^2 = gqdQ + gzdZ, 
respectively. To check that the system {dx,dy, ^1,^2} is not Frobenius, compute 
d^iA^iA^2=9Q{fuzfp-fpufz)dPAdQAduAdZ, 

_ ^ ^ _ , N ,n 7^ 7 iry mod dx,dy. (5.3) 

"6 A 4i A 42 = Jp[9uz9q - 9Qu9z)dP AdQ AduAdZ 

In each case, the coefficients on the right are nonzero, and we conclude that the transforma- 
tion is non-holonomic. 

5.3. Solving for Backlund transformations. In this subsection, we will set up systems of 
PDE whose solutions are Backlund transformations to the wave equation for some examples 
on the Goursat-Vessiot list. Although the existence of these transformations follows from 
the arguments of §4, here we will be able to go further in writing down explicit formulas 
for the transformations. Because we will work with specific Monge-Ampere equations on 
the list, we can take advantage of explicit formulas for the characteristic invariants. (These 
invariants are computed, for example, in the dissertation of M. Biesecker [I].) 
The general approach is as follows. We write a PDE on the list in the form 

s = F{x,y,u,p,q). (5.4) 

This form always has x and y as characteristic invariants, and we assume these are the 
only functionally independent invariants up to first order for the equation (i.e., we assume 
that the equation is not Monge-integrable). The Backlund transformation must take these 
invariants to corresponding characteristic invariants for the wave equation. By employing 
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a change of variables on the wave equation side, of the form X t— > Y h-> ipiY), and 

interchanging X and Y if necessary, we may assume that the transformation has 

x = X, y = Y. 

Now suppose that the remaining equations defining the Backlund transformation take the 
form 

P = fix,y,u,Z,P,Q), q = g{x,y,u,Z,P,Q). (5.5) 

The Monge-Ampere system on encoding the PDE fl5.4l) is generated algebraically by 
the contact form 6 = du — pdx — q dy and the 2-forms 

Qi = {dp — F{x, y, u, p, q)dy) A dx, VL2 = {dq — F{x, y, u, p, q)dx) A dy. 

The defining property of the Backlund transformation is that substituting fl5.5p into fii,il2 
must make them congruent to linear combinations of dP A dx and dQ A dy (the 2-forms 
defining the Monge-Ampere system for the wave equation) modulo 9 and the contact form 
on the wave equation side, 

9 = dZ — P dx — Q dy. 

In fact, Qi must become congruent to a multiple of dP A dx and Q2 congruent to a multiple 
of dQ A dy. Using (15.51) . we compute 

= {{fy + fug + fzQ - F)dy + fpdP + fQdQ) A dx mod 9, 9, 

and 

^2 = {{gx + Quf + 9zP - F)dx + gpdP + gqdQ) A dy mod 9, 9. 
We immediately conclude that fq = gp = 0, so that the transformation is of the form 

P = fix,y,u,Z,P), q = g{x,y,u,Z,Q), (5.6) 
and /, g must satisfy two additional first-order PDEs, 

fy = F{x,y,uJ,g) - fug - fzQ, (5.7) 
gx = F{x,y,uJ,g) - guf - gzP- (5.8) 

We derive additional first- and second-order PDEs that / and g must satisfy by differen- 
tiating the conditions so far. Taking derivatives with respect to Q in (15. 7p and P in (15.81) 
gives 

fz = {Fg-fu)gQ, 9z = iFp-gu)fp, (5.9) 
where the partials Fp = dF/ dp and Fg = dF/ dq are taken and then evaluated with p and q 
given by (15. 6p . As we will see in specific cases below, this will sometimes imply that / and 
g must be linear in P and Q. 

In what follows, let Ji and J2 denote the second-order characteristic invariants for the 
given PDE (whose existence makes the equation Darboux-integrable) , expressed in terms 
of x,y,u,p,q and the second-order jet coordinates r and t. (We make the convention that 
Ji is invariant along the characteristic curves where x is constant, and J2 is invariant when 
y is constant.) Then the Backlund transformation must take Ji and J2 to second-order 
characteristic invariants for the wave equation. In order to compute these additional con- 
straints, we must take total x- and ^-derivatives in (15. 6p to deduce how the second-order jet 
coordinates r and t transform in terms of those of the wave equation: 

r = fx + fuf + fzP + fpR, t = gy + gug + gzQ + 9qT. (5.10) 
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Requiring that, under these substitutions, Ji transforms to be a function of only x, P, R, 
and J2 transforms to be a function of only y, Q, T, will lead to additional second-order PDEs 
which / and g must satisfy. 

We now turn to specific examples. 

Equation IX (Liouville's equation). In this case, F = e", and the equations ( 15. 7p through 
(15. 9p become 

fy = e^ + UQgQ-g), (5.11) 

g. = e^ + 9u{Pfp-f), (5.12) 
fz = -fugq, (5.13) 
9z = -9ufp. (5.14) 

As mentioned in §1, the characteristic invariants are 

Ji=r- J2=t - \q^. 

Under (15.61) . the first invariant transforms as 

r - y = + fuf + fzP + fpR - \f. 

Requiring that this be a function of x, P, R only immediately implies that fp can depend on 
X and P only, and that the remaining terms have no dependence on m, Z or y. This gives us 
6 additional second-order PDEs for /: 

fpu = fpy = fpz = 0, d^, dy, dz{U + fuf + fzP - If) = 0. (5.15) 

Similarly, we also get 

gQu = QQy = Qqz = du,d:„dzigy + gu9 + 9zQ - =0- (5.16) 

Note that some of these second-order equations are redundant, in light of the derivatives of 
( KT3\f and flCTjl . 

Next, we derive additional equations by differentiation. Note that (15.111) shows that 
cannot be identically zero; then, taking a Q derivative of (15.131) shows that gqq = 0. Simi- 
larly, fpp = 0, so that / and g are linear in P and Q. Thus, we may set 

/(x, y, M, P, Z) = f^{x)P + f\x, y, u, Z), g{x, y, u, Q, Z) = g^{y)Q + g\x, y, u, Z). 

In particular, taking the terms in (15.151) . (I5.16P that are linear in P and Q respectively gives 

du {if - g')fl - ff) = 0, {{g' - f')gl - g'g') = 0. (5.17) 

Furthermore, equating the Z-derivative of (15.111) with the ^/-derivative of (I5.13p . and us- 
ing the M-derivatives of these equations to determine fyu and fzu, gives the compatibility 
condition 

i9'-f')f'u9l = 9X + e-g'; 

we similarly derive 

{f-9')fl9l = f'9l + e-f. 
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Adding and differentiating with respect to u, and using the values for and g^^ given by 
fl5.17p . shows that = —g^ = k for some nonzero constant k, and f^gi = Integrating 
the remaining equations shows that the most general form for the transformation is 

(5.18) 



, n ^ I'u + kZ + v(x) + w(y)\ , 
p = A;P + 2exp ( ^ ) +t;'(x), 



,^ i' u — kZ — v(x) — w(y)\ ,, , 
q = -kQ + exp ^ ^ - w'{y), 



where v{x),w{y) are arbitrary functions. 

Using the calculation (15.31) . it is easy to verify that none of these transformations is 
holonomic. 

In the next two examples, we will analyze the system of PDEs that / and g must satisfy 
using the techniques of exterior differential systems. 



Equation XIII. This PDE, 

s = F{u,x,y} := - — ■ — r^, 
[x + yY 

has second-order characteristic invariants 

Ji = r + — ^, J2 = t + 



X + y X + y 

in the x- and y-directions respectively (see [T], Appendix A). Substituting for p and r from 
(El]) and (15301) yields 

Ji= (f. + fuf + fzP+^)+fpR, 

so that fp and the expression in parentheses must be functions of x and P only. Similarly, 
we have 

J2 = (gy + gug + gzQ + — ^ ) + ggT, 

V x + yj 

hence gg and the expression in parentheses must be functions of y and Q only. 
In this case, (15.71) through (15. 9p specialize to 

fy = F-{g- QgQ)U, (5.19) 

g^ = F-{f- Pfp)gu, (5.20) 

fz = -ggfu, (5.21) 

gz = -fpgu. (5.22) 

If fu were identically zero, then fz would also be identically zero, but then fy = F = 
2u/{x + yY would give a contradiction. So, we may assume that and (similarly) g^ are 
nonzero on an open dense set. It then follows from (15.211) that gqq = and from (15.221) that 
fpp = 0, i.e., / and g are again linear in P and Q. 

Differentiating (I5.20I) . (I5.22I) with respect to x and Z, and equating mixed partials, enables 
us to solve for fp^ as 



2/ 



fp. = fuifp - gg) - "^^ (5.23) 

gu{x + yy 
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while from (I5.19p . fl5.2ip we similarly obtain 



= - M - (5.24) 

To encode the PDEs that / and g must satisfy as an exterior differential system, we will 
use x,y,u, Z, P,Q as independent variables, and use f, fx, fu, fp,ri and g, gy, gu, gQ,ti as 
dependent variables. (The role of the coefficients ri and ti will be made clear below.) We 
will regard these variables as coordinates on M}^. As stated above, we restrict to the open 
subset IX C M}^ where fp, fu, gq and gu are nonzero. 

The generator 1-forms are ipi through ip^, where 

ipi = -df + fxdx + fudu + fpdP + fzdZ + fydy, 
il)2 = -dg + gydy + gudu + gqdQ + gzdZ + gxdx, 
i'3 = -dfp + fpxdx, 
ipi = -dgq + gqydy, 

with fy, fz, gx, gz given by equations fl5.2ip through fl5.20p and fpx, gqy given by fl5.23p and 
fl5.24p . The remaining generators encode the rest of the condition that the second- 

order characteristic invariants be preserved. Differentiating the first term in Ji gives 

d (/. + + fuf + ^) - d{Q + Pdifz) + (d{fu) - j^^^ f 



Let 7]i be the 1-form on the right; then for any Backlund transformation rji must be a linear 
combination of dx and dP. In fact, since only the first term in rji can contain dP, the 
coefficient of dP in rji must be fpx- Thus, our remaining generators are 

i^5 = Vi- fpxdP - ridx, ip6 = r]2- gqydQ - hdy, 
where, based a similar calculation of dJ2, we set 

V2 = d{gy) + Qd{gz) + (^d{gu) - ^^^^^y ^ 9+ (^gu + (du^ - fpdJ + Fdx) . 

We seek to construct integral manifolds of the given differential ideal, i.e., submanifolds of 
U to which all the forms in the ideal pull back to be zero. An integral element for an EDS is 
an infinitesimal version of an integral manifold, i.e., a subspace in the tangent space to IX at 
some point, to which all the forms in the ideal restrict to be zero. Because we want integral 
manifolds which are graphs of functions of x, y, u, P, Q, Z, we will only consider integral 
elements which are 6-dimensional, and to which the differentials dx, dy, du, dP, dQ, dZ restrict 
to be linearly independent; we will call these admissible integral elements. 

Applying Cartan's Test to the Pfaffian system generated by ipi, ... ,ipQ shows that it has 
last nonzero Cartan character si = 4, but is not involutive, as the space of admissible inte- 
gral elements has 2-dimensional fiber at each point. However, the system becomes involutive 
after one prolongation, and this establishes the existence of the required Backlund transfor- 
mations. The last nonzero Cartan character of the involutive prolongation is Si = 2. By the 
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Cartan-Kahler Theorem (see [8], Chapter 7) we conclude that 6- dimensional integral sub- 
manifolds, satisfying the independence condition, exist through every point of IX, and that 
the construction of such submanifolds depends on a choice of 2 functions of one variable. 
The additional 1-forms that generate the prolongation include 

'/'s = (i(9j + yd + j dy + (/„9„ - J^J^'j t*^. 

which are actually defined on the original manifold M^^. These forms vanish on all integral 
elements of the original system, and if they had been included in the ideal, it would have 
been involutive with si = 2. 

The vanishing of ip-/, ips implies that fu and Qu are functions of x and y only. Moreover, 
forms ip7,ips define a smaller Pfaffian system, involving only fu,gu as functions of x and 
y, which satisfies the Frobenius condition. This means that fu{x,y) and guix,y) can be 
determined by solving systems of ODE. Once these are determined, substituting the solutions 
into f l5.23p and f l5.24p gives a Frobenius system which may be solved for the functions fp{x) 
and ggiy). Then / and g may be determined by integrating first-order PDE, with / including 
an arbitrary function of x and g an arbitrary function of y. 

For example, by observing that + gu must satisfy a Riccati equation as a function of 
X + y, we are led to a solution 



{x + y)x' " {x + y)y' 

Substituting these into fl5.23p . fl5.24p leads to fp + gq = k{x + y)/{xy) for a constant k. It is 
simplest to choose k = with fp = l and gq = —1. Integrating then gives the solution 

x{x + y) y[x + y) 

Proposition 5.1. All Backlund transformations between (XIII) and the wave equation are 
holonomic. 

Proof. As noted in §5.2[ the holonomic condition is equivalent to the Pfaffian system on 
spanned by dx, dy and 

dp - fudu = {-gqfu)dZ + fpdP, dq - g^du = {-fpgu)dZ + gqdQ, mod dx, dy 

being Frobenius. It is straightforward to check that d{fpdP—gqfudZ) and d{gqdQ — f pg^dZ) 
are zero modulo dx, dy and the 1-forms of the above EDS. (For example, d{fp) = modulo 
dx, dy, ipi, ■ ■ ■ , '08) and the same is true for d{gq) and d{fu)-) □ 



Equation IV. This PDE has the form 

u 

where a and (5 are arbitrary solutions of the ODE df/dt ±t/f = K for some fixed 7^ 0. 
(We will take the plus sign in these equations, the computation for the other sign being 
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completely analogous.) In this case, (15. 9p takes the form 



fz={^(K--^)-^]gQ, (5.25) 



u V Pig) 

9z=[^(K-47^]-9u]fp. (5.26) 



u \ a{f) 

(From now on, instead of writing a{f) and P{g), a and /3 will be understood to be composed 
with /(x, y, u, Z, P) and g{x, y, u, Z, Q) respectively.) 

Unlike in previous examples, here it is not valid to conclude that / and g are linear in P 
and Q. In fact, differentiating (15.251) with respect to Q gives 

= {u^U - [KP - g)a)gQQ + ^ (5.27) 

enabling us to determine gqQ. (If the coefficient in front were identically zero, then f3{q) 
would be identically equal to a constant times q, which contradicts K ^ 0.) Similarly, 
differentiating (I5.26P yields 

= (nag. - (Ka - f)P)fpp + - - f)f) f2_ (5.28) 

The characteristic invariants for (IV) are 

r a t (3 

Ji = , = "5 • 

a u p u 

Substituting for p and r from (15. 6p and (l5.1Up gives 

J, = ^±M±M±M_«, (6.29) 

a u 

so that fp/a must be a function of x and P only. Setting the derivatives of this with respect 
to M, and Z equal to zero yields 

/^ fpfu J. ( f\ fpfy J. f f\ fpfz 



fPu=[K-^^, fpy=(K~^^, fpz=[K-^ (5.30) 
\ a / a \ a J a \ a J a 

where fy is given by (15.70 and fz is given by (15.250 . Similarly, from the T coefficient in J2 
we get that gq/P must be a function of y and Q only, and hence 

f 9\ 9q9u ( 9\ 9q9x [ 9\ 9q9z 

We may also differentiate (15.251) and (15.261) to obtain equations for fp^ and gqy. 

We encode the various first- and second-order partial differential equations for / and g 
derived so far into an exterior differential system generated by 1-forms ipi, . . . ,1^^, as we did 
for equation (XIII). Unlike the previous example, we do not need to prolong, but instead 
obtain integrability conditions which take the form 

«5'n = Pfu, (5.32) 

and 

{upfu - {.K(3 - g)a)gQ = {uag^ - {Ka - f)(3)fp. (5.33) 
(Note that, by using (I5.25p . (I5.26p . this implies that agz = Pfz ) With these conditions 

incorporated into the EDS, it becomes involutive with last nonzero character si = 2. 
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Solutions of this system may be obtained by observing that the quantities 



A = — , /i = — , 7 = a — (A/ + nP)u, S = P — {\g + nQ)u (5.34) 
a a 

must be functions of x, y, u and Z only, and satisfy the following compatible system of 
first-order PDE: 



{K-u\)X, 



dX 

du 


= -uX^ + KX^ - -A, 
u 




du 
du 


{u^ 


A2 - KuX + 1) 

u 


1^7 
7 du 


ld5 
5 du 


dX 


d^ 




(9/i 


= {K- 


uX)fi'^, 


1 ^7 


1 d5 


dZ 


du'' 




'dZ 


7aZ 


~ 6dZ 


dX 


{u^x^ - Kux + 1; 


1 

-7, 


d^ 


= {K- 


uX) — , 
u 


dS 


= {K- 


dx 


u^ 


dx 


dx 


dX 


{u^x^ - Kux + r 




(9/i 


= {K- 


mA)— , 
u 


d^ 


d5 


dy 




dy 


dy 


dx 



u 



A solution (7, 6, A, /i) to this PDE system may be constructed by integrating successively 
in the w-direction, the Z-direction, the x-direction and the ?/-direction. (Note that the x- 
dependence of 7 and the ^/-dependence of 6 are given by arbitrary functions.) Once 7 and 6 
are known, they implicitly determine / and g. 

Using (15. 3p . one can check that the resulting Backlund transformations are holonomic if 
and only if, in the above system, m^A^ — KuX -|- 1 = 0. Thus, holonomic transformations 
exist, and depend on fewer arbitrary constants but the same number of arbitrary functions. 
For example, if K = 2, then a solution to the above system is given by 

1 —1 v'(x)u ^ w'(y)u 



u Z — v{x) — w{yy Z — v{x)—w{y)' Z — v{x) ~ w{y) 

Then, using (15.341) . a holonomic Backlund transformation is implicitly defined by 

{v'jx) - P)u {w'{y) - Q)u 

a{p)~p = - — — , f3(^q)-q- 



Z — v{x) — w{y) Z — v{x) — w{y) 



5.4. Summary. Besides equations (IV) and (XIII) discussed above, we have also investi- 
gated the exterior differential system for Backlund transformations to the wave equation for 
equations (V), (VII), (IX), (XI) and (XII). Even if explicit formulas are not available, in each 
case we use the Cartan-Kahler Theorem to determine (in terms of the last nonzero Cartan 
character) the size of the solution set, in both the holonomic and non-holonomic cases. The 
results are summarized in the table below. 
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Equation 


Monge-Integrable 


Holonomic BTs 


Non-holonomic BTs 


I 


no 


yes, si 


= 2 


yes. 


Sl = 


2 


II 


no 


no 




yes. 


Sl = 


2 


III 


no 


no 




yes. 


Sl = 


2 


IV 


no 


yes, si 


= 2 


yes. 


Sl = 


2 


V 


no 


yes, si 


= 2 


yes. 


Sl = 


2 


VI 


yes 


yes, si 


= 3 


yes. 


S3 = 


1 


VII 


no 


no 




yes. 


Sl = 


2 


VIII 


yes 


yes, si 


= 3 


yes. 


S3 = 


1 


IX 


no 


no 




yes. 


Sl = 


2 


X 


yes 


yes, si 


= 3 


yes, S3 = 


1 


XI 


no 


yes, si 


= 2 




no 




XII 


no 


yes, si 


= 2 




no 




XIII 


no 


yes, si 


= 2 




no 





Note that the approach described in §5.3l is not feasible for the Monge-integrable equations 
(VI, VIII and X), but for completeness we include them in the table, together with the results 
from the analysis in §4.21 The Cartan character for the system for holonomic Backlund 
transformations for such equations is variable, depending on whether one considers the cases 
(1), (2), or (3), as described in U3. 

It is interesting to note that equations (I), (IV) and (V) have both holonomic and non- 
holonomic transformations, in roughly the same degree of generality. In fact, it is possible 
that these two kinds of Backlund transformations linking the same pair of equations may be 
closely related. In our previous paper [3], we pointed out that the transformation (Z.I) is 
a composition of two simpler transformations, a holonomic Backlund transformation to the 
wave equation, and a contact transformation from the wave equation to itself. It is possible 
that, more generally, the non-holonomic transformations for these equations are obtainable 
from holonomic transformations in this way. 

6. Concluding Remarks 

In this section, we will indicate some interesting directions in which the results in this 
paper might be extended, and some important questions about Backlund transformations to 
which the techniques in this paper may be relevant. 

(1) The set of equivalence classes (under contact transformations) of second-order Monge- 
Ampere equations to which the results of §3 in this paper apply is relatively small, 
confined to the equations on the Goursat-Vessiot list. It would be interesting to see if 
the arguments in that section could be applied to hyperbolic systems of class A; > 1. 
In other words, given a hyperbolic system I of class k, linked to the standard wave 
equation by a Backlund transformation, can one prove that the prolongation of I is 
Darboux-integrable? Likewise, given a hyperbolic EDS I of class k, such that its pro- 
longation is Darboux-integrable, does there exist a Backlund transformation between 
X and the Monge-Ampere system for the standard wave equation? (The argument 
given at the end of Chapter 7 in [S] shows that there is a Backlund transformation 
between the wave equation and the prolongation of X; however, for practical purposes 
it is desirable to have a Backlund transformation between systems of as low an order 
as possible, so that one has a smaller system of ODE to solve in order to construct 
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solutions.) These hyperbolic systems would include, for example, the Monge- Ampere 
equations which are Darboux-integrable at third order, which have not been classified 
and are thought to comprise a much larger set. 

(2) It is a theorem of Sophus Lie that no Monge- Ampere equation of the form u^y = f{u) 
is Darboux-integrable (after arbitrary many prolongations) except when f{u) = 
exp{au + h) for constants a and h (see [6], Chapter IX). Consequently, important 
equations like sine-Gordon cannot have a Backlund transformation to the wave equa- 
tion. Instead, the Backlund transformation (11.11) for sine-Gordon produces solutions 
to the same PDE as we started with. This is known as an auto- Backlund transfor- 
mation^ such transformations play an important role in the theory of completely 
integrable PDE [9]. 

It is therefore of interest to try to identify those Monge- Ampere equations which 
have non-trivial auto-Backlund transformations. We remark that for such transfor- 
mations, the Monge- Ampere systems X on M and X on M must be contact-equivalent, 
i.e.,. there must be a diffeomorphism $ : M — * M which pulls back X to X. 

vr / \ 7f 

/ $ \_ 
M 

Necessary conditions for the existence of such a diffeomorphism may be derived from 
the fact that it is required to preserve the differential invariants of the Monge- Ampere 
systems. (See 0, §2.1, for a derivation of these invariants using the method of 
equivalence.) 

(3) Our previous paper began the exploration of parametric Backlund transforma- 
tions using the method of equivalence. Such transformations contain an arbitrary 
parameter in the Backlund system; for example, an arbitrary nonzero parameter A 
may be interpolated in the sine-Gordon auto-Backlund transformation (11.11) to give 

Vx-u^ = ^ sin((M + v)/2), 

^y + Uy = sin((M - v)/2). 

One observes that this system differs from (11.11) merely by scaling x by A and y by 
A^^ — a change of variables which is a symmetry of the sine-Gordon equation but 
not of the system (11.11) . This scaling symmetry can also be applied to the Backlund 
transformation (11.21) . to produce a parametric transformation 

Zx = Ux- 2Aexp((u + z)/2), Zy = -Uy + ^exp((M - z)/2), 

where u{x,y) satisfies Liouville's equation and z{x,y) solves the wave equation. (In 
fact, this transformation is derived from the most general form (15.181) by setting 
k = 1 and choosing v{x) = 2 In A and w{y) = 0.) In |4j it is shown that these 
transformations can be generated from a non-parametric Backlund transformation 

^This terminology is not universally accepted; Hongyou Wu [12] has proposed that transformations be- 
tween different PDEs be known as Miura transformations, and the term Backlund transformation be reserved 
for what we are calling auto-Backlund transformations. 
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by starting with a symmetry vector field on M, choosing a hft into which is 
not a symmetry of the Pfaffian system 3, but such that pulhng 3 back to S x M 
via the 1-parameter family of diffeomorphisms generated by the lift gives a family 
of transformations. The same approach can be taken with other transformations 
discussed in §4.1; for example, the transformation (IZ.Iip may be generalized to a 
parametric transformation 



ZZr — XuUx = \/ Z^ — Am^a/A + z\^ ZZy — XuUy = V z"^ — ^ \ + z^ 



for A > 0, where u satisfies (II) and z solves the wave equation. (This is obtained by 
starting with the symmetry of (II) that simultaneously scales m, x and y.) 

With these examples in evidence, and given the importance of parametric Backlund 
transformations in the study of 'soliton' equations, it is desirable to characterize those 
transformations that may be made to depend on an arbitrary parameter by lifting 
symmetry vector fields. 

References 

[I] M. Biesecker, Geometric Studies in Hyperbolic Systems in the Plane, Ph.D. thesis, Utah State 
University, 2003. 

[2] R. Bryant, P. Griffiths, D. Grossmann, Exterior differential systems and Euler-Lagrange partial 

differential equations. University of Chicago Press, 2003. 
[3] J. CleUand, Homogeneous Backlund transformations of hyperbolic Monge- Ampere systems, 

Asian J. Math. 6 (2002), 433-480. 
[4] J. CleUand, T. Ivey, Parametric Backlund transformations I: Phenomenology, Trans. Amer. 

Math. Soc. 357 (2005), f06f-f093. 
[5] R. Gardner, The Method of Equivalence and Its Applications, CBMS-NSF Regional Conf. Ser. 

in Appl. Math. 58, SIAM, Philadelphia, 1989. 
[6] E. Goursat, Logons sur I'integration des equations aux derivees partielles du second ordre, 

Gauthier-Villars, 1896-1898. 
[7] -, Recherches sur quelques equations aux derivees partielles du second ordre, Annates de la 

Faculte de Toulouse (deuxieme serie) 1 (1899), 31-78. 
[8] T. Ivey, J.M. Landsberg, Cartan for Beginners: Differential geometry via moving frames and 

exterior differential systems (with J.M. Landsberg), Graduate Studies in Mathematics vol. 61, Amer- 
ican Mathematical Society, 2003. 
[9] A. P. Fordy, A Historical Introduction to Solitons and Backlund Transformations, in "Harmonic Maps 

and Integrable Systems" (ed. Fordy and Wood), Vieweg, 1994. 
[10] M. Juras, Geometric Aspects of Second-Order Partial Differential Equations in the Plane, 

Ph.D. thesis, Utah State University, 1997. 

[II] E. Vessiot, Sur les equations aux derivees partielle du second ordre integables par la methode 
de Darboux, J. Math Pures Appl. 18 (1939), 1-61 and 21 (1942) 1-66. 

[12] Xifang Cao, Hongyou Wu, Chuanyou Xu, On Miura transformations among nonlinear partial 
differential equations, J. Math. Phys. 47 (2006), no. 8. 

[13] M.Yu. Zvyagin, Classification of Backlund transformation of Second-Order Partial Differen- 
tial Equations, Mat. Zametki 29 (1981), 829-842 (EngHsh translation: Math. Notes 29 (1981), 422-429). 

[14] -, Second order equations reducible to z^y = by a Backlund transformation, Dokl. Akad. 
Nauk SSSR 316 (1991), 36-40 (English translation: Soviet Math. Dokl. 43 (1991) 30-34). 

Dept. of Mathematics, 395 UCB, University of Colorado, Boulder, CO 80309-0395 
E-mail address: Jeanne.Clelland@colorado.edu 

Dept. of Mathematics, College of Charleston, 66 George St., Charleston SC 29424-0001 
E-mail address: IveyT@cofc.edu 



48 



